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ABSTRACT 

Aims. A significant fraction of early-type galaxies (ETGs) show emission lines in their optical spectra. We aim at understanding the 
powering mechanism and the origin of the ionized gas in ETGs, and its connection with the host galaxy evolution. 
Methods. We analyzed intermediate-resolution optical spectra of a sample of 65 ETGs, mostly located in low density environments 
and biased toward the presence of ISM traces, for which we already derived in the previous papers of the series the stellar population 
properties. To extract the emission lines from the galaxy spectra, we developed a new fitting procedure that properly subtracts the 
underlying stellar continuum, and that accounts for the uncertainties due to the age-metallicity degeneracy. The emission line lumi- 
nosities derived in annuli of increasing galacto-centric distance were used to constrain the excitation mechanism and the metallicity 
of the ionized gas. 

Results. Optical emission lines are detected in 89% of the sample. The detection fraction drops to 57% if only the galaxies with 
EWQia + [NII])>3 A are considered. The incidence and strength of emission do not correlate either with the E/SO classification, 
or with the fast/slow rotator classification. Comparing the nuclear i<rjl6 emission with the classical [0III]/Hj8 vs [NII]/Hq' di- 
agnostic diagram, the galaxy activity is so classified: 72% of the galaxies with emission are LINERs, 9% are Seyferts, 12% are 
Composite/Transition objects, and 7% are non-classified. Seyferts have young luminosity-weighted ages (< 5 Gyr), and appear, on 
average, significantly younger than LINERs and Composites. Excluding the Seyfert galaxies from our sample, we find that the spread 
in the ([OIII], Ha or [Nil]) emission strength increases with the galaxy central velocity dispersion cr^, low-cr^ galaxies having all 
weak emission lines, and high-cr^ galaxies displaying both weak and strong emission lines. Furthermore, the [NII]/Hq- ratio tends 
to increase with cr^. A spatial analysis of the emission line properties within the individual galaxies reveals that the [NII]/Hq' ratio 
decreases with increasing galacto-centric distance, indicating either a decrease of the nebular metallicity, or a progressive "softening" 
of the ionizing spectrum. The average nebular oxygen abundance is slightly less than solar A comparison with the stellar metallicities 
derived in Paper III shows that the gas oxygen abundance is » 0.2 dex lower than that of stars. 

Conclusions. The stronger nuclear (r< r^/lG) emission can be explained with photoionization by PAGE stars alone only for a; 22% of 
the LINERs/Composite sample. On the other hand, we can not exclude an important role of PAGE star photoionization at larger radii. 
For the major fraction of the sample, the nuclear emission is consistent with excitation from a low-accretion rate AGN, fast shocks 
(200 -500 km/s) in a relatively gas poor environment (n < 100 cm"^^) , or coexistence of the two. The derived [SII]67 17/6731 ratios 
are consistent with the low gas densities required by the shock models. The derived nebular metallicities suggest either an external 
origin of the gas, or an overestimate of the oxygen yields by SN models. 

Key words. Galaxies: elliptical and lenticular, cD - Galaxies: ISM - Galaxies: active-Galaxies:abundances. 



1. Introduction 

Early-type galaxies (ETGs) have long been considered to be in- 
ert stellar systems, essentially devoid of gas and dust. However, 
this view has radically changed since a number of imaging and 
spectroscopy studies from both the ground and space have re- 

* Based on observations obtained at the European Southern 
Observatory, La Silla, Chile. 

** Full Tables 2 and 3 are only available at the CDS via anonymous 
ftp to cdsar c . u- strasbg . f r 



vealed the presence of a multiphase interstellar medium (ISM): 
a hot (T ~ 10*' - 10^ K), X-ray emitting halo (F orman & 
Jones [T9851 Fabbiano et al. IT9921 O' Sullivan et al. l200TT i: a 
warm (T ~ 10'^ K) gas (often referred to as "ionized gas") 
(e.g. Philli ps et a l. 11^861 Sadler & Gerhard 1985; see also 
Goudfrooij 119991 for a review); and colder components de- 
tected in the mid and far infrared, CO and HI (see review 
of Sadler et al. l2002l Temi et al. l2004l Br essan et al. l2006l 
Sage & Welch 2006; Morganti et al l2006l Helmboldt [20071 
Oosterloo et al. ,2007. Panuzzo et al. I2007I I. The hot gas com- 
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ponent dominates the ISM in E/SO galaxies, and only compara- 
tively small quantities are detected in the warm and cold phases 
(Bregman et al. 1992). The derived masses of the warm ISM 
range from ~ 10^ to - W^Mq (Phillips et al. 1986; Trinchieri 
& di Serego Alighi eri |199ll: Macchetto et al. ,1996. Gou dfrooij 
et al. |1994a|[T9Pblll994ct Goudfrooij & de Jong [T995] l. orders 
of magnitude below the X-ray emitting hot ISM and the stellar 
mass. Ionized gas is detected in 40-80% of early-type galaxies 
via its optical emission lines (Caldwell 1984; Phillips et al. 119861 
Macchetto et al.lT996t Sarzi et al. i2006. ,2009; Yan et al. l2006l 
Serra et al. 120081) . The ionized gas is also usually morphologi- 
cally associated with dust (Kim 1989', Goudfrooij et al. 1994c'). 

Despite the number of studies, several issues remain still 
open. The first question is the origin of the ISM in ETGs. Narrow 
band imaging centered around Ha'+[NII]/16584 shows that the 
ionized gas presents a variety of morphologies, from regular, 
disk-like structures, to filamentary structures (Demoulin-Ulrich 
et al. [T984; Buson et al. [1993. Macchetto et al.il99 6; Zeilinger 
et al. [19961 Martel et al. [20041 Sarzi et al. l2006b . Evidence for 
the acquisition of external gas comes from kinematical studies 
showing that the angular momentum of the gas is often decou- 
pled from that of the stars (Bertola et al. 119921 van Dokkum 
& Franx il995 ; Caon et al. 2000) . even if according to Sarzi et 
al. ( 120061 ) the angular momenta are inconsistent with a purely 
external origin for the gas. The fact that the ionized gas emission 
is always associated with dust (e.g. Tran et al. 12001]) tends to 
exclude "cooling flows" as the origin for the warm gas in E/SO 
galaxies (Goudfrooij 1 1 9991 ) . In the scenario proposed by Sparks 
et al. ( 119891 ) and de Jong et al. ( 119901 ). the dusty filaments can 
arise from the interaction of a small gas rich galaxy with the gi- 
ant elliptical, or from a tidal accretion event in which the gas 
and dust are stripped from a passing spiral (see e.g. Domingue 
et al.laop. Tal et al. [20091) . 

The history of star formation and evolution leaves its chemi- 
cal signature in the ISM of ETGs. Historically, X-ray abundance 
measurements of the hot ISM have been problematic, as typified 
by the so-called Fe discrepancy (Arimoto et al. 119971 ). More re- 
cently, Chandra and XMM measurements of X-ray bright galax- 
ies at the center of groups have supported roughly solar or even 
slightly supersolar abundances for the ISM (e.g. Gastaldello 
& Molendi i2002t Buote et al. [20031 Tamura et al. |2003|), and 
have claimed no evidence for very subsolar Fe abundances 
(Humphrey & Buote [20061 Ji et al. 120091 ). These studies also re- 
port anomalously low oxygen abundances compared to iron and 
magnesium. This is in conflict with results from SN stellar yields 
if one assumes that the source of ISM is the material injected into 
interstellar space by evolved stars. Up to now, the only study of 
metallicity in the warm ISM is that of Athey & Bregman (2009 ), 
who derived oxygen abundances for seven ETGs from optical 
emission lines. The authors find an average solar metallicity, fa- 
voring an internal origin of the warm ISM. However, this study 
is limited by the small sample. Determining the metallicity of the 
warm ISM is fundamental to create a link between the hot gas 
phase and the galaxy stellar population, and to discriminate be- 
tween an internal and an external origin of the interstellar matter 
in E/SO galaxies. 

The second still open issue concerns the ionizing source of 
the warm gas. Optical spectroscopic studies show that ETGs 
are typically classified as LINERs (Low-Ionization Nuclear 
Emission-line Regions, Heckman I 19801 ) according to their emis- 
sion line ratios (e.g. Phillips et al. '1986; Yan et al. i2006l) . 
However, there is still strong debate about the ionization 
mechanism in LINERs. At present, the most viable excitation 
mechanisms are: a low accretion-rate AGN (e.g.. Ho il999bt 



Kewley et al. 120061 Ho I2009al ). photoionization by old post- 
asymptotic giant branch (PAGB) stars (e.g. Trinchieri & di 
Serego AHghieri 1991 ; Binette et al. [19941 Stasins ka et al . 20081, 
fast shocks (Koski & Osterbrock 1976; Heckman [19801 Dopita 
& Sutherland 1995. Allen et al. ,2008) . The AGN mechanism is 
strongly supported by the detection of broad emission lines in 
optical specti-a (e.g. Ho et al. ll997at Wang & Wei 2008 ), by the 
properties of the X-ray emission (Terashima et al. 20(30; Flohic 
et al. 120061 Gonzalez-Martin et al. 120091), by the presence of UV 
and X-ray variability (e.g Maoz et al. 120051 Plan et al. 120101), 
and by the evidence that massive black holes (BHs) appear to 
be a generic component of galaxies with a bulge (Magorrian 
et al. [19981 Ho ll999al Kormendy 2004). PAGB stars are good 
candidates because their radiation, much harder than that from 
young stars, is able to reproduce the observed emission-line ra- 
tios in LINERs (Trinchieri & di Serego Alighieri ll991t Binette 
et al. |1994t Stasinska et al. 120081) . The most compelling ev- 
idence in support of the PAGB mechanism is the finding that 
the emission-line flux correlates very well with the host galaxy 
stellar luminosity within the emission-line region (Macchetto et 
al. 1996), and that the emission-line flux distribution follows that 
of the stellar continuum (Sarzi et al. 2006, 2009 ). These corre- 
lations suggest that the source of ionizing photons is distributed 
in the same way as the stellar population. LINERs emission has 
also been observed in extranuclear regions associated with large- 
scale outflows and related shocks (Lipari et al. 2004), or regions 
shocked by radio jets (Cecil et al. 2000 ). 

Evidence is growing that more than one mechanism may 
be at work in LINERs. In a recent study, Sarzi et al. ( 120091 ) 
showed that the role of AGN activity is confined to the central re- 
gions, whereas at larger radii the stellar and nebular fluxes follow 
each other, thus suggesting a stellar ionizing source. Eracleous 
et al. (120101 ) showed that low accretion-rate AGNs do not pro- 
duce enough ionizing photons to explain the observed Ha lumi- 
nosities within a 2" x 4" aperture, so that other mechanisms are 
needed. 

With the aim of understanding the formation history of 
ETGs, and the connection between the stellar populations and 
the ionized gas, we have analyzed intermediate-resolution opti- 
cal spectra for a sample of 65 field E/SO galaxies (Rampazzo 
et al. 12005! hereafter Paper I; Annibali et al. ,2006. hereafter 
Paper II). The stellar population properties (i.e. ages, metallic- 
ities, and elements abundance ratios) were derived at different 
radii in Annibali et al. (2007) (hereafter Paper III) through the 
analysis of the Lick indices. In this paper we present the study of 
the optical emission lines in annuli of increasing galacto-centric 
distance. The purpose is to characterize the properties of the ion- 
ized gas, its origin, and the possible excitation mechanisms at 
different radii. 

The paper is organized as it follows. In section 2 we pro- 
vide a brief overview of the sample. In section 3 we describe the 
method used to extract the emission-line fluxes from the optical 
spectra. In section 4 we provide a spectral classification through 
standard diagnostic diagrams, and investigate possible correla- 
tions with the host galaxy properties. In section 5 we determine 
the oxygen metallicity. In section 6 we compare our data with the 
models. A summary and the conclusions are given in section 7. 

2. Sample Overview 

The Rampazzo et al. ( 2005!) + Annibali et al. (^2006) sample 
(hereafter R05+A06 sample) was selected from a compilation of 
ETGs showing ISM traces in at least one of the following bands: 
IRAS 100 yum. X-ray, radio, HI and CO (Roberts et al. , 19911 . 



Table 1. Sample overview 
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Galaxy Id 


RSA 


RC3 


o-c 


re 




Pxyz 


Age 


Z 


[a/Fe] 








km s"' 


arcsec 


km s"' 


gal Mpc~' 


Gyr 






NGC 128 


S02(8) pec 


SO pec sp 


183 


17.3 


4227 




9.7 ± 1.7 


0.024 ± 0.004 


0.16 ± 0.03 


NGC 777 


El 


El 


317 


34.4 


5040 




5.4 ±2.1 


0.045 ± 0.020 


0.28 ± 0.10 


NGC 1052 


E3/S0 


E4 


215 


33.7 


1475 


0.49 


14.5 ± 4.2 


0.032 ± 0.007 


0.34 ± 0.05 


NGC 1209 


E6 


E6: 


240 


18.5 


2619 


0.13 


4.8 ± 0.9 


0.051 ± 0.012 


0.14 ± 0.02 


NGC 1297 


S02/3(0) 


SABO pec: 


115 


28.4 


1550 


0.71 


15.5 ± 1.2 


0.012 ± 0.001 


0.29 ± 0.04 


NGC 1366 


E7/S01(7) 


SOsp 


120 


10.6 


1310 


0.16 


5.9 ± 1. 


0.024 ± 0.004 


0.08 ± 0.03 


NGC 1380 


S03(7)/Sa 


SAO 


240 


20.3 


1844 


1.54 


4.4 ± 0.7 


0.038 ± 0.006 


0.24 ± 0.02 


NGC 1389 


S01(5)/SB01 


SAB(s)0-: 


139 


15.0 


986 


1.50 


4.5 ± 0.6 


0.032 ± 0.005 


0.08 ± 0.02 


NGC 1407 


E0/S01(0) 


EO 


286 


70.3 


1766 


0.42 


8.8 ± 1.5 


0.033 ± 0.005 


0.32 ± 0.03 


NGC 1426 


E4 


E4 


162 


25.0 


1443 


0.66 


9.0 ± 2.5 


0.024 ± 0.005 


0.07 ± 0.05 


NGC 1453 


EO 


E2 


289 


25.0 


3906 




9.4 ±2.1 


0.033 ± 0.007 


0.22 ± 0.03 


NGC 1521 


E3 


E3 


235 


25.5 


4165 




3.2 ± 0.4 


0.037 ± 0.006 


0.09 ± 0.02 


NGC 1533 


SB02(2)/SBa 


SBO- 


174 


30.0 


773 


0.89 


11.9 ±6.9 


0.023 ± 0.020 


0.21 ±0.10 


NGC 1553 


S01/2(5)pec 


SA(r)0 


CO 


65.6 


1280 


0.97 


4.8 ± 0.7 


0.031 ± 0.004 


0.10 ± 0.02 


NGC 1947 


S03(0) pec 


SO- pec 


142 


32.1 


1100 


0.24 


5.9 ± 0.8 


0.023 ± 0.003 


0.05 ± 0.02 


NGC 2749 


E3 


E3 


248 


33.7 


4180 




10.8 ± 2.3 


0.027 ± 0.006 


0.25 ± 0.04 


NGC 2911 


S0porS03(2) 


SA(s)0: pec 


235 


50.9 


3131 




5.7 ± 2.0 


0.034 ± 0.019 


0.25 ±0.10 


NGC 2962 


RSB02/Sa 


RSAB(rs)0-H 




23.3 


2117 


0.15 








NGC 2974 


E4 


E4 


220 


24.4 


1890 


0.26 


13.9 ± 3.6 


0.021 ± 0.005 


0.23 ± 0.06 


NGC 3136 


E4 


E: 


230 


36.9 


1731 


0.11 


1.5 ±0.1 


0.089 ± 0.004 


0.36 ± 0.02 


NGC 3258 


El 


El 


271 


30.0 


2778 


0.72 


4.5 ± 0.8 


0.047 ± 0.013 


0.21 ± 0.03 


NGC 3268 


E2 


E2 


227 


36.1 


2818 


0.69 


9.8 ± 1.7 


0.023 ± 0.004 


0.34 ± 0.04 


NGC 3489 


S03/Sa 


SAB(rs)+ 


129 


20.3 


693 


0.39 


1.7 ±0.1 


0.034 ± 0.004 


0.05 ± 0.02 


NGC 3557 


E3 


E3 


265 


30.0 


3038 


0.28 


5.8 ± 0.8 


0.034 ± 0.004 


0.17 ± 0.02 


NGC 3607 


S03(3) 


SA(s)0: 


220 


43.4 


934 


0.34 


3.1 ± 0.5 


0.047 ± 0.012 


0.24 ± 0.03 


NGC 3818 


E5 


E5 


191 


22.2 


1701 


0.20 


8.8 ± 1.2 


0.024 ± 0.003 


0.25 ± 0.03 


NGC 3962 


El 


El 


225 


35.2 


1822 


0.32 


10.0 ± 1.2 


0.024 ± 0.003 


0.22 ± 0.03 


NGC 4374 


El 


El 


282 


50.9 


1060 


3.99 


9.8 ± 3.4 


0.025 ± 0.010 


0.24 ± 0.08 


NGC 4552 


SOl(O) 


E 


264 


29.3 


322 


2.97 


6.0 ± 1.4 


0.043 ± 0.012 


0.21 ± 0.03 


NGC 4636 


E0/S01(6) 


EO-1 


209 


88.5 


937 


1.33 


13.5 ±3.6 


0.023 ± 0.006 


0.29 ± 0.06 


NGC 4696 


(E3) 


E-i-1 pec 


254 


85.0 


2958 


0.00 


16.0 ± 4.5 


0.014 ± 0.004 


0.30 ±0.10 


NGC 4697 


E6 


E6 


174 


72.0 


1241 


0.60 


10.0 ± 1.4 


0.016 ± 0.002 


0.14 ± 0.04 


NGC 5011 


E2 


El-2 


249 


23.8 


3104 


0.27 


7.2 ± 1.9 


0.025 ± 0.008 


0.25 ± 0.06 


NGC 5044 


EO 


EO 


239 


82.3 


2704 


0.38 


14.2 ± 10. 


0.015 ± 0.022 


0.34 ±0.17 


NGC 5077 


S01/2(4) 


E3+ 


260 


22.8 


2764 


0.23 


15.0 ± 4.6 


0.024 ± 0.007 


0.18 ± 0.06 


NGC 5090 


E2 


E2 


269 


62.4 


3421 




10.0 ± 1.7 


0.028 ± 0.005 


0.26 ± 0.04 


NGC 5193 


SOl(O) 


Epec 


209 


26.7 


3711 




6.8 ± 1.1 


0.018 ± 0.002 


0.26 ± 0.04 


NGC 5266 


S03(5) pec 


SAO-: 


199 


76.7 


3074 


0.35 


7.4 ± 1.4 


0.019 ± 0.003 


0.15 ± 0.05 


NGC 5328 


E4 


El: 


303 


22.2 


4671 




12.4 ± 3.7 


0.027 ± 0.006 


0.15 ± 0.05 


NGC 5363 


[S03(5)] 


10: 


199 


36.1 


1138 


0.28 


12.1 ±2.3 


0.020 ± 0.004 


0.16 ± 0.05 


NGC 5638 


El 


El 


168 


28.0 


1676 


0.79 


9.1 ± 2.3 


0.024 ± 0.008 


0.24 ± 0.05 


NGC 5812 


EO 


EO 


200 


25.5 


1930 


0.19 


8.5 ±2.1 


0.027 ± 0.008 


0.22 ± 0.05 


NGC 5813 


El 


El-2 


239 


57.2 


1972 


0.88 


11.7 ± 1.6 


0.018 ± 0.002 


0.26 ± 0.04 


NGC 5831 


E4 


E3 


164 


25.5 


1656 


0.83 


8.8 ± 3.5 


0.016 ±0.011 


0.21 ± 0.09 


NGC 5846 


SOl(O) 


E0+ 


250 


62.7 


1709 


0.84 


8.4 ± 1.3 


0.033 ± 0.005 


0.25 ± 0.03 


NGC 5898 


S02/3(0) 


EO 


220 


22.2 


2267 


0.23 


7.7 ± 1.3 


0.030 ± 0.004 


0.10 ± 0.03 


NGC 6721 


El 


E-I-: 


262 


21.7 


4416 




5.0 ± 0.8 


0.040 ± 0.007 


0.24 ± 0.02 


NGC 6758 


E2 (merger) 


E-I-: 


242 


20.3 


3404 




16.0 ± 2.5 


0.016 ± 0.002 


0.32 ± 0.05 


NGC 6776 


El pec 


E+pec 


242 


17.7 


5480 




2.7 ± 0.5 


0.033 ± 0.010 


0.21 ± 0.05 


NGC 6868 


E3/S02/3(3) 


E2 


277 


33.7 


2854 


0.47 


9.2 ± 1.8 


0.033 ± 0.006 


0.19 ± 0.03 


NGC 6875 


S0/a(merger) 


SAB(s)0- pec: 




11.7 


3121 










NGC 6876 


E3 


E3 


230 


43.0 


3836 




9.8 ± 1.6 


0.023 ± 0.003 


0.26 ± 0.03 


NGC 6958 


R?S01(3) 


E+ 


223 


19.8 


2652 


0.12 


3.0 ± 0.3 


0.038 ± 0.006 


0.20 ± 0.03 


NGC 7007 


S02/3/a 


SA0-: 


125 


15.4 


2954 


0.14 


3.4 ± 0.6 


0.031 ± 0.010 


0.15 ± 0.05 


NGC 7079 


SBa 


SB(s)0 


155 


19.8 


2670 


0.19 


6.7 ± 1.1 


0.016 ± 0.003 


0.21 ± 0.05 


NGC 7097 


E4 


E5 


224 


18.1 


2404 


0.26 


10.5 ± 2.4 


0.024 ± 0.005 


0.30 ± 0.05 


NGC 7135 


SOI pec 


SAO- pec 


231 


31.4 


2718 


0.32 


2.2 ± 0.4 


0.047 ± 0.010 


0.46 ± 0.04 


NGC 7192 


S02(0) 


E-I-: 


257 


28.6 


2904 


0.28 


5.7 ± 2.0 


0.039 ± 0.015 


0.09 ± 0.05 


NGC 7332 


S02/3(8) 


SO pec sp 


136 


14.7 


1207 


0.12 


3.7 ± 0.4 


0.019 ± 0.002 


0.10 ± 0.03 


NGC 7377 


S02/3/Sapec 


SA(s)0-i- 


145 


36.9 


3291 




4.8 ± 0.6 


0.020 ± 0.002 


0.10 ± 0.03 


IC 1459 


E4 


E 


311 


34.4 


1659 


0.28 


8.0 ± 2.2 


0.042 ± 0.009 


0.25 ± 0.04 


IC 2006 


El 


E 


122 


28.6 


1350 


0.12 


8.1 ± 0.9 


0.026 ± 0.003 


0.12 ± 0.02 


IC 3370 


E2 pec 


E2+ 


202 


38.6 


2934 


0.20 


5.6 ± 0.9 


0.022 ± 0.004 


0.17 ±0.04 


IC 4296 


EO 


E 


340 


41.4 


3762 




5.2 ± 1.0 


0.044 ± 0.008 


0.25 ± 0.02 


IC 5063 


S03(3)pec/Sa 


SA(s)0+: 


160 


26.7 


3402 
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IC 1459 (H) 



[on] [Nelll] H/S [OIII][OIII] [NII]Ha[NII] [SII][SII 




3700 3800 3900 4900 5000 6500 6600 6700 



m 

NGC 5813 (W) 




3700 3800 3900 4900 5000 6500 6600 6700 



A (A) 

Fig. 1. Central (r <r(,/16) spectra for a high emission galaxy (IC 1459) and for a weak emission galaxy (NGC 5813) in the wavelength 
ranges where the most important emission lines are measured. For each galaxy, we show in the top panels the observed spectrum 
(black line), the underlying stellar spectrum model (red), and the fit to the observed spectrum obtained by adding emission lines 
(Gaussian curves) to the underlying stellar continuum (blue). In the bottom panels we show the subtracted spectrum, where the 
emission lines emerge. Notice that in this figure the subtracted stellar continuum is an SSP of age, Z, and extinction equal to the 
average of all the solution with;^f^ < 2. However, this is only for illustration purposes, and it is slightly different from the procedure 
described in Section 3.1. 



All galaxies belong to the Revised Shapley Ames Catalog of 
Bright Galaxies (RSA) (Sandage & Tammann ll987l) and have 
a redshift lower than ~ 5500 km s Because of the selection 
criteria, the sample is biased toward the presence of emission 
lines. Table 1 summarizes the main characteristics of the sample. 
Column (1) gives the galaxy identification name; Col. (2) and 



(3) provide the galaxy morphological classification according to 
RSA (Sandage & Tammann 119871 1 and RC3 (de Vaucouleurs 
et al. 119911), respectively: only in few cases do the two cata- 
logues disagree in the distinction between E and SO classes; 
Col. (4) gives the central velocity dispersion from HYPERLEDA 
(http://leda.univ-lyonl.fr//) ; Col. (5) gives the effective radius 
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from RC3; Col. (6) gives the galaxy systemic velocity, 
which is lower than ~5000 km s"' ; Col. (7) provides the richness 
parameter p vyj (Tully 1988 ): it represents the density of galaxies 
brighter than -16 apparent B-mag in the vicinity of the entry, in 
galaxies xMpc"-'. Cols (8), (9) and (10) Hst the age, metallicity 
and [a/Fe] ratios derived in Paper III. 

The galaxies of the R05-HA06 sample are mainly located 
in low density environments. The local density of the galaxies 
varies from ~ 0.1 Mpc"^, typical of very isolated galax- 
ies, to pivj a; 4 Mpc^^*, characteristic of denser regions in the 
Virgo cluster. For comparison, in the Tully ( 1988 ) catalogue ob- 
jects like NGC 1399 and NGC 1389, members of the Fornax 
cluster, have pj:v;=l-59 and 1.50 Mpc"-', respectively. Thus, the 
R05-HA06 sample, even though biased towards low density en- 
vironments, contains a fraction of galaxies in relatively dense 
environments. The sample spans a wide range in central velocity 
dispersion, a; from 115 to 340 km s"' (see Papers I and II for 
details). 

Following the RC3 classification, the sample is composed of 
~ 70% ellipticals and ~ 30% lenticulars. However, if we classify 
the galaxies according to their amount of rotation (see Appendix 
A), we end up with ~ 70% fast rotators (F) and ~ 30% slow 
rotators (S). 

Intermediate resolution (FWHM a!7.6 A at 5550 A) spectra 
in the (3700 - 7250) A wavelength range were acquired with the 

I. 5 m ESO-La Silla telescope for the 65 galaxies of the sam- 
ple. The spectra were extracted for 7 apertures of increasing ra- 
dius (1.5", 2.5", 10", r,/10, r,/8, r,/4, and r,/2), corrected for 
the galaxy ellipticity, and 4 adjacent annulii (r <r(,/16, r(,/16 < r 
<Xej%, rf/8< r <r(,/4, and rp/4 < r <Xel'T)- The data reduction and 
the computation of the Lick indices are described in Papers I and 

II. In Paper III we derived ages, metallicities, and [cc/Fe] ratios 
by comparing the data with our new Simple Stellar Population 
(SSP) models. No stellar population parameters are derived for 
IC 5063, since line emission is too strong to allow a measure- 
ment of the Lick indices. Summarizing the results of Paper III, 
we derive a large age spread, with SSP-equivalent ages ranging 
from a few Gyrs to a Hubble time. The galaxies have metal- 
licities and [c/Fe] ratios above solar. Both the total metallicity 
and the [a/Fe] present a positive correlation with the central ve- 
locity dispersion, indicating that the chemical enrichment was 
more efficient and the duration of the star formation shorter in 
more massive galaxies. We also find that the youngest objects 
in our sample are all located in the lowest density environments 
(piyz < 0.4 Mpc~^). We suggest that the young galaxies in the 
lowest density environments underwent secondary episodes of 
star formation, which we call "rejuvenation episodes". Within 
the individual galaxies, the stellar metallicity tends to decrease 
from the center outwards. We derive an average metallicity gra- 
dient of AlogZ/Alog(r/r,) ~ -0.21. 

3. Emission lines 

3. 7. Starlight Subtraction and Line Intensities 

There are a number of emission lines that we expect to detect 
within the wavelength range (3700 - 7250) A sampled by our 
spectra: [Oil] ^3727, [Nelll] /13869, H/3 /14861, [OIII] AAA959, 
5007, [OI] /16300, Ha A 6563, [Nil] ^^6548, 6584, [SII] AA6111, 
6731. 

The emission line fluxes were measured in residual spec- 
tra obtained after subtracting the stellar population contribution 
from the observed galaxy spectra. The starlight subtraction is a 
crucial step for the correct determination of the emission line 




-1x10-1= lO-is 2x10-15 

Flux H/3(erg cm"^) 



Fig. 2. Ha versus li/3 emission fluxes for the galaxy sample mea- 
sured in a central r <ijl6 region. IC 5063, with derived emis- 
sion fluxes falling outside the plot boundaries, is not included. 
Full dots denote the seven "no emission" galaxies (NGC 5638, 
NGC 5831, NGC 1407, NGC 1366, NGC 1426, NGC 5812, and 
NGC 1389) in the sample. The solid lines indicate regions of 
different reddening. 



properties, in particular for HjS and Ha, which are superimposed 
to the Balmer absorption features of the underlying stellar pop- 
ulation. The H/3 emission, much fainter than the Ha iFH/3 ~ 
1/3 FHa in absence of extinction), is the most difficult to mea- 
sure. With only few exceptions, we never see in our spectra a 
true emission feature in Hy6, but rather an infilling of the stel- 
lar absorption line. It follows that a reliable measurement of this 
line rests on the appropriate modeling of the underlying starlight 
contribution. 

To subtract the underlying stellar population from the galaxy 
spectra we used new SSPs (Bressan, unpublished, see also 
Clemens et al. |2009l Chavez et al. [20091 ) . The SSPs were com- 
puted from the isochrones of Padova 94 (Bertelli et al. 119941 ) 
with the revision of Bressan, Granato & Silva (1998), including 
a new AGB mass-loss treatment calibrated on LMC cluster col- 
ors. These SSPs are particularly suited for our analysis because 
they make use of the MILES spectral library (Sanchez-Blazquez 
et al. 12006'). MILES consists of 2.3 A FWHM optical spectra of 
=» 1000 stars spanning a large range in atmospheric parameters, 
and represents a substantial improvement over previous libraries 
used in population synthesis models. 

For our study, we considered the galaxy spectra extracted 
in 4 concentric annulii (r <r<,/16, ie/l6 < r <r(,/8, rj,/8< r <r(,/4, 
and Te/4- < r <rJ2), as described in Section 2. The galaxy con- 
tinua were fitted with the new SSPs, previously smoothed to 
match the instrumental resolution and the galaxy velocity dis- 
persion, through a criterion. In the fit, we considered se- 
lected spectral regions chosen to be particularly sensitive either 
to age or to metallicity. For each SSP of given age and metallic- 
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For each line, we adopted the average emission flux com- 
puted as: 



Fig. 3. E(B - V) distribution derived for our sample from the 
pHa/pH^ ratios. The shaded histogram denotes SO galaxies, 
while the dotted line is for E galaxies. Vertical arrows indicate 
the median values. 



ity, and E(B - V) ranging from the Galactic foreground value to 
E(B-V)max = 0.5,the;r' is: 



1 



Ni- ^ (Ni)i 

^ ~ Nr ^ (N^)i ^ 



(FobsWj) - FsSpC/lj) X Cnorm.i)^ 



(1) 



where Nr is the number of selected spectral regions, (N/l)i 
is the number of points in the /-th region, Fobs(/lj) is the ob- 
served galaxy flux at a given wavelength, Fssp(/lj) x Cnorm.i is 
the SSP flux normalized to the observed flux in the /-th region, 
cr(/lj) is the error on the observed flux. The Nr spectral regions 
have widths between ^ 100 A and ^ 300 A. The normalization 
constant Cnorm.i is computed in a 20 A wide region adjacent 
to the /-th band. This approach guarantees that the solution is 
mostly driven by the strengths of the absorption features, rather 
than by the global slope of the observed galaxy spectrum. This 
minimizes the effect of extinction and of the uncertainties due 
to flux calibration on the solution. The spectral regions were se- 
lected in order to include both features particularly sensitive to 
age (the relative strength of the Ca II H -H K lines, the 4000 A 
break (D4000) , die Hy and H(5 lines), and features more sensitive 
to metaUicity (the Fe lines at A 4383, 4531, 5270, 5335, and the 
Mg absorption features around A 5175). The region around Hy6, 
where emission is expected, was obviously excluded from the fit. 
On the other hand, negligible emission is expected in the higher- 
order Balmer lines (Hy, and He, blended with the Call H line 
at our resolution). The residual spectrum around each line is de- 
rived by normalizing the SSP in two continuum bands adjacent 
to the line of interest. 

Because of the degeneracy between age, metaUicity, and ex- 
tinction, fits with similarly good values can produce signifi- 
cantly different residual spectra, in particular around HyS. Thus, 
to obtain a statistically meaningful determination of the emission 
lines, we considered all the N fits with;^f^ < 2, and computed the 
emission lines on the corresponding residual spectra. 



y2<2 



XF; 



Z,t,Av 



x'-<2 



(2) 



where Fz.t,Av is line flux obtained by subtracting the SSP of 
metaUicity Z, age t, and reddening Ay to the observed galaxy 
spectrum, and e"-*^' is its weight. The line flux Fz,t,Av is deter- 
mined by fitting the residual spectrum with Gaussian curves of 
variable width and intensity. The emission line width is treated 
as a free parameter because the velocity dispersion of the gas 
can be significantly different from that of the stars. More specif- 
ically, single Gaussians were used to fit all the lines, with the 
exceptions of [SII]/l/l6717, 31, which was fitted with a sum of 2 
Gaussians, and the [NII]/1/16548, 84 + Ha complex, which was 
fitted with a sum of 3 Gaussians. The [OII]/13727 feature is actu- 
ally a doublet (/1/13726, 29), but the two lines appear completely 
blended in our spectra. 

The error on the average emission flux is: 



X-<'2 



x(F: 



Z,t,Av 



(3) 



This error takes into account the uncertainty due to the in- 
trinsic degeneracy of the stellar populations. Because of the high 
signal-to-noise of our spectra, this is actually the largest source 
of error in the measure of the emission fluxes. 

The use of a single SSP in the fitting procedure may lead to 
spurious effects where recent star formation is present in some 
amounts. We showed in Paper III that some galaxies in our sam- 
ple have very young luminosity weighted ages (<4 Gyr), possi- 
bly as a consequence of "rejuvenation" episodes. Therefore, we 
repeated the fitting procedure using, instead of one SSP, a com- 
bination of a young (t < 8 Gyr) and an old (t> 8 Gyr) SSP. This 
approach guarantees a more reliable spectral subtraction under 
the Balmer lines for the "rejuvenated" galaxies. The two-SSP fit- 
ting procedure results in some differences for the galaxies with 
faint Hy6 emission. This may be important when discriminating 
between LINERs and Seyfert through the [OIII]/Hy6 ratio (see 
Section 4). On the other hand, we checked that there are no sig- 
nificant differences in the resulting emission lines for galaxies 
randomly selected among those without signs of recent star for- 
mation. 

In Fig. [U we show the stellar continuum subtraction pro- 
cedure for two illustrative examples: IC 1459, with relatively 
strong emission lines, and NGC 5813, with weak emission lines. 
Notice that in the figure the subtracted stellar continuum is an 
SSP of age, metaUicity, and extinction equal to the average of 
all the N solution with;^^^ < 2. This is slightly different from our 
procedure, where we have subtracted the individual solutions N 
times from the same galaxy spectrum, and have computed the 
final emissions as average of the N emission fluxes. In Fig. [T] 
the emission features are easily visible in IC 1459. The errors 
on the emission line fluxes introduced by the uncertainties in 
the stellar continuum subtraction are thus relatively small. In the 
case of NGC 5813, instead, the determination of the emissions 
is dominated by starlight subtraction effects, in particular for HjS 
and Ha. Consequently the errors on the emission lines are quite 
large. 

We provide in Table 2 the emission line fluxes in units of 
10"'^ erg s"' cm"^ arcsec"^ derived for the R05-i-A06 sample in 
4 annuli of increasing galacto-centric distance (r <Xell6, re/16 < 
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Fig. 4. li/3, [Olll]A 5007, Ha and [Nil] A 6584 EW distributions at r <rJ16 for the total galaxy sample, for the SO and E subsamples, 
and for the fast (F) and slow (S) rotators subsamples. IC 5063, with EW(I^) ~ 15, EW([OIII]) ~ 122, EWCHo-) ~ 58, and EW([NII]) 
~ 38, falls outside the plot boundaries. The EWs are not corrected for extinction. The vertical arrows indicate the median values for 
the E, SO, F, and S subsamples. 




r <re/8, Te/^< r <re/4, and 1^/4- < r <j:J2). In Table 3 we provide 
the line equivalent widths (EW), derived as F^/Fcont, where Fcom 
is the flux computed in a continuum band adjacent to the line of 
interest. The high errors in HyS are due to its intrinsic faintness, 
accompanied by a strong dependence of the H/? absorption line 
on age. The values given in Tables 2 and 3 are not corrected 
for reddening. The extinction corrected emissions can be easily 
derived from the E(B - V) values given in Col. 3 of Table 2, and 
computed as described in the next paragraph. 

3.2. Extinction 

The extinction in the lines was derived through the relative ob- 
served strengths of the Balmer lines. The intrinsic flux ratio 
(FHa/FH/?),) is * 2. 85 for HII regions and ^ 3.1 for AGN-like 
objects (Osterbrock 1 1 9891 1. Values larger than these are due to 
absorption from an intervening medium. 

More specifically, given an observed FHa/FH/? ratio, the 
reddening is computed as: 



E(B - V) = 



logio [(FHa/FH;8)/(FHg/FHyg)o] 
0.4 X Rv X (Aj^g/Av - Ah„/Av) 



(4) 



where A^ is the magnitude attenuation at a given wave- 
length, and A^g/Av, Aua/Ay, and Ry = Av/E(B - V), depend 
on the adopted extinction curve. We adopted the Cardelli et 
al. (1989) extinction law, where Ry = 3.05. We also assumed 
(FHa/FHy6)o = 3.1, since all the galaxies in the R05+A06 sam- 
ple display AGN - like emission-line ratios (see Section 4). The 



assumption of (FHa/FHy3)o = 2.85 does not change significantly 
our results. 

In Figure |2] we plot the Ho- versus H/? fluxes derived in 
a central r <Te/l6 region. The majority of the galaxies have 
E(B - V) < 0.3. Some galaxies display instead very strong ex- 
tinction, as high as E(B - V) a; 1 .5 or even more. They have very 
faint H/3 intensities, accompanied by relatively strong Ho- emis- 
sions. Vice versa, a few objects fall below the E(B - V) = line, 
in the region of "negative reddening". In these cases, the emis- 
sion derived in H/3 is incompatible (too large) with the derived 
Ha. An obvious explanation is that the underlying continuum 
was not subtracted properly. More specifically, the fact that we 
underestimate the Ha/H/3 ratio suggests that we have subtracted 
populations that are systematically too young or too metal poor 
To verify this hypothesis, we simulated the situation in which 
a too-young population is subtracted from the galaxy spectrum 
in the following way. We created a synthetic spectrum by adding 
H/3 and Ha emission lines (in the ratio FHa = 3.1 FH/3) to an old 
SSR Then we re-derived the emission intensities by subtracting 
from the synthetic spectrum increasingly younger SSPs. We ob- 
tained that the FHa/FH/3 ratio decreases as the A age increases. 
The effect is as much stronger as the initial emission intensities 
are weaker It is due to a combination of the stronger dependence 
of the H/3 absorption line on age, and of the intrinsic faintness of 
the H/3 emission line with respect to Ha. 

Finally, there are some galaxies whose HyS and Ho- emissions 
are both consistent with zero within a 3 cr error (NGC 1366, 
NGC 1389, NGC 1407, NGC 1426, NGC 5638, NGC 5812, and 
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Fig. 5. Emission EW derived in annuli of increasing galacto- 
centric distance (r <rf/16, re/16 < r <re/8, re/8< r <re/4, and 
re/4 < r <re/2) over the central value (r <re/16) as a function 
of r/rg. The top and bottom panels refer to the [Nil] 6584 A and 
Ha lines, respectively. The shaded regions indicates the average 
trends and the + 1 o" levels. Only galaxies classified as H or W in 
Table 4 were included. Galaxies with flat or increasing emission 
trends were labeled. 



NGC 5831). These galaxies show no emission in the other lines 
as well, and are classified as no emission (N) galaxies. 

The derived E(B - V) values at different galacto-centric dis- 
tances are given in Col. 3 of Table 2. For the galaxies with 
"negative reddening", and for the galaxies with E(B - V) lower 
than the Galaxy foreground extinction, the E(B - V) was set to 
the foreground value. Very large E(B - V) values are derived 
for NGC 777, NGC 1521, NGC 3136, NGC 3557, NGC 6776, 
NGC 7007, NGC 3489. Such sti'ong extinction values, not ob- 
served in the continuum, suggest that the distribution of the dust 
is patchy. Notice however that the errors in the derived E(B - V) 
values are quite large, due to the high uncertainty in the H/3 emis- 



sion. The reddening distributions for the E and SO classes are 
shown in Fig. [3] The median E(B - V) values for the two classes 
are 0.11 and 0.26 respectively. However, a Mann-Withney U test 
shows that there is not a significant difference between the two 
distributions. 

Once computed the E(B - V) values, we corrected the ob- 
served emission fluxes: 



Fobs X 10" 



where 

A, 



A = |^|x3.05xE(B-V) 



(5) 



(6) 



is taken from the extinction curve, and A is the wave- 
length of the line of interest. 

3.3. Results 

The galaxies are classified in Table 4 according to the inten- 
sity of their nuclear emission lines. We identify three cate- 
gories: no emission galaxies (N), where no emission lines are 
detected within a 3 cr error; weak emission line galaxies (W), 
with EW(Hq' -I- [NII]6584)<3 A; and strong emission line galax- 
ies (H), with EWCHa + [NII]6584)>3 A. 

Emission lines are detected in 58 out of 65 galaxies (~ 89% 
of the sample), while strong emission is present in 57% of the 
sample. In the central r <re/16 region, we derive median EWs of 
-0.51, 0.72, 1.4 and 1.8 for H/3, [OIII]/15007, Ha, and [Nll]A 
6584, respectively. For the same lines, the median percentage 
errors amount toa;31%,ll%,6%, and , 1 %. 

We investigated possible differences in the emission inten- 
sity between the E and SO subsamples, and between the fast (F) 
and slow (S) rotators subsamples. The distributions and the rela- 
tive median values are shown in FigH] An appreciable difference 
is observed only between F and S for the median [Nil] value. 
However, Mann- Whitney U tests provide P > 0.05 (for a two- 
tailed test) for all the distributions, indicating that there is no 
significant difference between the E, SO, F and S subsamples. 

To quantify the importance of the emission with respect to 
the underlying stellar continuum at increasing galacto-centric 
distance, we studied the behavior of the [NII]/16584 and of the 
Ha lines in the 4 annuli defined in Section 2. Notice that the 
use of adjacent annuli instead of "apertures" allows us to isolate 
the true emission contribution from different regions within the 
galaxy. The [NII]/16584 line is relatively stronger, and thus more 
easily measured also in the external annuli. Furthermore, it is 
also very poorly affected by uncertainties in the underlying stel- 
lar population subtraction, and it has the lowest percentage error. 
On the other hand the Ha presents the advantage that, being a re- 
combination line, directly reflects the ionizing photon flux, but 
its measure is much more affected than the [Nil] by uncertainties 
in the underlying population subtraction. The behaviors of the 
observed [NII]/([NII]<re/16) and Ha/(Ha<re/16) ratios with in- 
creasing galacto-centric distances are shown in Fig.|5] Negative 
values are consistent with no emission. The plotted values were 
not corrected for reddening, but extinction could affect the EWs 
if the nebular emission was more reddened than the continuum. 
Indeed, Tran et al. (120011) showed that the presence of emission 
lines in ETGs tends to be associated with the presence of clumpy 
dust. 

Fig. |5] shows that, with some exceptions (NGC 3489, 
NGC 5898, NGC 7007, NGC 7079, NGC 7192, NGC 7377), 
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Fig. 6. Extinction corrected BPT diagnostic diagram for the galaxy sample in annuli of increasing galacto-centric distance (r <rf/16, 
ie/l6 < I <rf/8, r(,/8< r <ie/4, and 1^/4 < r <iJ2). Triangles are for galaxies with very large reddenings (E(B - V) > 1.5 at r <rf/16). 
The solid curve is the "maximum starburst line" of Kewley et al. (120011) . while the dashed line indicates the empirical division 
between pure star-forming galaxies and AGN-HII composite (or transition) objects from Kauffmann et al. (120031 1. The horizontal 
line at log([OIII]/Hj6) = 0.5 separates "Seyfert" and LINERs galaxies (Kewley et al. l2006l l. 



the emission EWs tend to decrease from the center outwards. 
The average decrease in the [Nil] EW is ^ 0.7, 0.5, and 0.4 of 
the central value at ~0.1 r^,, ~0.2 r^, and ~0.4 r^,, while the Ha 
EW decreases more gently (~ 0.8, 0.7, and 0.6 at the same radii). 
The steeper decrease of the [Nil] compared to the Ha indicates 
a progressive softening of the ionizing spectrum from the galaxy 
center to more external regions. The decrease in the Ho- EW im- 
plies that the ionizing flux decrease more rapidly than the stellar 
continuum around ~ 6500 A. However, given the uncertainties, 
we can consider significant only the decrease from the center to 
~ 0.1 I'f,. Deeper observations are needed to trace with higher 
confidence the EW behavior out to larger radii. 



The galaxies NGC 7079, NGC 3489, NGC 7007, NGC 7 192, 
NGC 7377, NGC 5898 deviate from the average decreasing 
trend, presenting flat or increasing emission EWs (although the 
Ha of NGC 3489, NGC 5898 and NGC 7079 drops in the most 
external annulus). NGC 3489 and NGC 7007 have large red- 
denings, thus their increasing trend could be due to stronger ex- 
tinction toward the galaxy center. Unfortunately, the errors in 
the derived E(B - V) values are large, and it is difficult to es- 
tablish from our data if indeed the reddening decreases with 
galacto-centric distance. On the other hand, the other galaxies 
(NGC 7079, NGC 7192, NGC 7377, and NGC 5898) show very 
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Fig. 7. Distribution of emission line ratios over the central r 
<re/16 value for the galaxies in the four BPT diagrams of Fig. |6] 
The dashed, empty, and full histograms correspond respectively 
to the rf/16 < r <vJ8, tJSk r <tJ4, and rJA < r <tJ2) annuli. 
Vertical arrows indicate the average values. From the galaxy cen- 
ter outwards, the average [N1I]/Hq' / [Nlll/Hac ratios drop to 
0.85, 0.73, and 0.69 in the three annuli, respectively. For the 
same annuli, we derive average [OIII]/Hy8 / [OIII]/HjSc values 
of 0.92, 0.77, and 0.64. 

little or almost null extinction, thus the observed trends are likely 
to be intrinsic. 

4. Spectral classification 

4.1. Diagnostic Diagrams 

Galaxies were classified through the standard [Olll]/I5007/ li/3 
versus [NII]/16584/ Ha diagnostic diagram (Baldwin, Phillips 
and Terlevich (119811 hereafter BPT). Using a combination of 
stellar population synthesis models and photoionization mod- 
els, Kewley et al. ( 1200 11 1 identified a "maximum starburst line" 
in the BPT diagram. Galaxies lying above this hne are likely 
to be dominated by an AGN. Kauffmann et al. (I2003I I later re- 
vised the Kewley et al. (2001) classification scheme introduc- 
ing an empirical line to divide pure star-forming galaxies from 
AGN-Hll Composite objects (also called "Transition objects"), 
with significant contribution from both AGN and star formation 
(see also Stasiriska et al. 120061 ). The optical spectra of composite 
objects can be due to either 1) a combination of star formation 
and a Seyfert nucleus, or 2) a combination of star formation and 
LINERs emission (Kewley et al. l2006b . 

From the distribution in the BPT diagram of the galaxies 
from the Sloan Digital Sky Survey, Kewley et al. ( I2006I I in- 
troduced an empirical horizontal line to separate "Seyfert" and 
LINERs. Even though LINERs show [Nil]/ Ha ratios larger than 
the star-forming galaxies, indicating a harder ionization contin- 
uum, it is still under debate if they are powered by active galactic 
nuclei or not (see Section 1). 

The BPT diagram for the R05-t-A06 sample is shown in 
Fig. |6] The [Nll]/}6584/ Ha line ratio was preferred to the 



Table 4. Classification from optical lines in the nuclear (r <re/16) 
region 



Galaxy Id 


Emission Type" 


Activity Class_^ 




TJ 

ri 




NGC 777 


W 


Sy/ LINER 


NGC 1052 


H 


LINER 


NGC 1209 


H 


LINER 


NGC 1297 


H 


LINER 


NGC 1366 


N 




NGC 1380 


H 


LINER 


NGC 1389 


N 




NGC 1407 


N 




NGC 1426 


N 






rl 


T Twm? 

i^iiN 1j,K 


NGC 1521 


w 


LINER 


NGC 1533 


H 


LINER 


NGC 1553 


W 


LINER 


NGC 1947 


H 


LINER 


NGC 2749 


H 


LINER 


NGC 2911 


H 


LINER 


NGC 2962 


H 




NGC 2974 


H 


LINER 


NGC 3136 


H 


LINER/ Sy 


NGC iZjo 


H 


Comp 


NGC 3268 


H 


LINER 


NGC 3489 


H 


Sy/ LINER 


NGC 3557 


W 


LINER 


NGC 3607 


H 


LINER 


NGC 3818 


Traces 




NGC 3962 


H 


LINER 


NGC 4374 


H 


LINER 


NGC 4552 


W 


Comp 


NGC 4636 


H 


LINER 


NGC 4090 


H 


LINER 


NGC 4697 


W 


LINER 


NGC 5011 


W 


LINER 


NGC 5044 


H 


LINER 


NGC 5077 


H 


LINER 


NGC 5090 


H 


LINER 


NGC 5193 


W 


Comp 


NGC 5266 


H 


LINER 


NGC 5328 


W/Traces 


Comp 


NGC 5363 


H 


LINER 


NGC jd3o 


N 




NGC 5812 


N 




NGC 5813 


W 


LINER 


NGC 5831 


N 




NGC 5846 


H 


LINER 


NGC 5898 


W 


LINER 


NGC 6721 


W/Traces 


Comp 


NGC 6758 


H 


LINER 


NGC 6776 


H 


LINER/ Sy 


NGC 6868 


H 


LINER 


NGC 05/ J 


W 




NGC 6876 


W /Traces 


Comp 


NGC 6958 


H 


Sy/ LINER 


NGC 7007 


W 


Sy/ LINER 


NGC 7079 


W 


LINER 


NGC 7097 


H 


LINER 


NGC 7135 


H 


LINER 


NGC 7192 


W 


LINER 


NGC 7332 


W/Traces 




NGC 7377 


W 


LINER 


IC 1459 


H 


LINER 


IC 2006 


W 


Comp 


IC 3370 


H 


LINER 


IC 4296 


H 


LINER 


IC 5063 


H 


Sy 



" N=no emission; W=weak emission (EW(H(!r + [NII]6584) <3 A); 
H=high emission (EW(Hq' -i- [NII]6584) >3 A) 
* At larger radii, some LINERs are classified as Composites. 
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log([NII]/Ha) 

Fig. 8. BPT diagrams of some galaxies representing the general trends of the whole sample in the four annuli: r <ic/16 (red cross), 
Te/i€> < r <rf/8 (blue open circle), r(./8< r <re/4 (green full square), and iJ4 < r <re/2 (yellow full triangle). Notice that while the 
error in [OIII] / li/3 is large, the error in [NII]/Ha is significantly smaller. 



[SII](/16717 + /i6731)/ Ha and [OI]^6300/ Ha ratios because 
of the larger signal to noise of the [Nil] line in our spectra. 
Emission line ratios were corrected for extinction adopting the 
E(B - V) values derived in Section 3.2 and reported in Col. 3 of 
Table 2, even though the BPT diagram is almost insensitive to 
reddening. The galaxies with negative reddening, or reddening 
lower than the foreground value, where corrected adopting the 
galactic foreground reddening from NED . 

According to their nuclear (r < re/16) emission, the major- 
ity of our galaxies are classified as LINERs. A few galaxies 
(NGC 3258, NGC 4552, NGC 5193, NGC 5328, NGC 6721, 
NGC 6876, IC 2006) fall in the region of "Composites", and 
possibly contain a combined contribution from both star for- 
mation and AGN. IC 5063, NGC 777, NGC 3489, NGC 7007, 
and NGC 6958, fall in the Seyfert region. Among them, only 
IC 5063 can be classified as a bona fide Seyfert galaxy. The oth- 
ers all have very high uncertainties in the [OIII]/ H^S ratios, due 
to the low H/3 emission (nuclear EW<0.3 A), and are consistent 
with a LINERs classification within the errors. These galaxies 
also display very large reddenings (E(B - V) >1). Vice versa, 
NGC 3136 and NGC 6776 are LINERs, but are consistent with 
being Seyferts within the errors. Among the Seyferts, a spectral 
classification is present in the literature only for NGC 3489 (Ho 
et al.' 1997b!, Sarzi et al. 2006), NGC 777 (Ho et al. 1997b), and 
NGC 6958 (Saraiva et al. 2001). These studies classify the first 
two galaxies as Seyferts, and the last one as LINER. 



The final spectral classification obtained from the nuclear (r 
<iJ16) lines is given in Col. 3 of Table 4. 

4.2. Gradients 

We investigated how the distribution of the galaxies in the BPT 
diagram changes at annuli of increasing galacto centric distance 
(r <r,V16, rjl6 < r <iJ8, iJ8< r <rJ4, and r,V4 < r <r,V2). This 
emerges by comparing the four panels of Fig. |6l Only a few 
galaxies have sufficiently high S/N ratios that allows to measure 
the emission line ratios in the most external annuli. 

The clear result is that, from the center outwards, the bulk of 
the galaxies moves left-down in the BPT diagram. An increas- 
ing number of galaxies pass from LINERs to "Composite" ob- 
jects as the emission lines are measured at larger galacto-centric 
distances. Likely, this effect would not have been observed by 
considering apertures instead of annuli, since the central emis- 
sion is almost always dominating. The eflrect is better quanti- 
fied in Fig. |7] where we plot the distribution of the [NII]/Hq' 
and [OIII]/H/3 ratios over the central r <rf/16 values (C) for the 
different annuli. The peak of the ([NII]/Ha)/([NII]/Ha)c distri- 
bution moves from 1 to values < 1 at larger radii. The aver- 
age ([NII]/Ha')/([NII]/Ha)c values are ^ 0.85, 0.73, and 0.69 at 
rp/16 < r <i'e/^, i'e/^< r <rp/4, and re/4 < r <i'e/2, respectively. 
These correspond to horizontal shifts of ^ -0.07, -0.14, and 
-0.16 dex in the BPT diagram. However, given the large er- 
rors in the emission line measures in the most external annu- 



12 



Annibali et al.: Nearby early-type galaxies with ionized gas. 




Fig. 9. Ha and [Nil] equivalent widths, [NII]/Hq' ratio, and [OIII]/HyS ratio versus galaxy central velocity dispersion 0-^, age, metal- 
licity Z, and [a/Fe] ratio. In the two top rows, full dots are for fast rotators (F), while empty dots are for slow rotators (S). In the 
two bottom rows, full dots are for high emission galaxies (H), and empty dots are for weak emission galaxies (W), as defined in 
Section 3 and Table 4. For the third row panels, we show least square fits for the H subsample. 



' ' ' ' I ' ' ' ' I ' ' ' ' I ' 1' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' '1 ' I ' ' ' I ' ' ' I ' ' ' 




0.6 



Fig. 10. Age, metallicity, and [a/Fe] distributions for Seyfert-like galaxies (S), LINERs (L), and "Composites" (C). 
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log([NII]/Ha) 

Fig. 11. BPT diagnostic diagram for Seyfert galaxies (left panel), "Composites" (center), and LINERs (right panel). The data were 
color-coded according to the galaxy age as derived in Paper III: blue stars are for ages younger than 4 Gyr, green squares for ages 
between 4 and 8 Gyr, and red circles for ages older than 8 Gyr 



lus, we consider this result significant only out to r <rf/4. The 
([OIII]/Hy8)/([OIII]/Hy8)c distribution is less peaked, in part also 
because of the large errors in the derived [OIII]/Hy6 ratios. The 
average values are 0.92, 0.77, and 0.64. for the three annuli, cor- 
responding to vertical shifts of =s -0.04, -0.11, and -0.19 dex 
in the BPT diagram. However, given the intrinsic faintness and 
large errors of the [OIII] and HyS lines compared to [Nil] and Ha, 
none of our conclusions will be based on the spatial behavior of 
the [OIII]/H/3 ratio. 

When the galaxies are inspected individually, we see that 
for the majority of them the [NIIJ/Hof ratio decreases mono- 
tonically with galacto-centric distance. The only exceptions 
are: IC 5063, NGC 1521, NGC 3136, NGC 3489, NGC 5813, 
NGC 5898, NGC 6721, NGC 6876, NGC 7135 and NGC 7192. 
In all these galaxies the [NIIl/Ha ratio does not have a clear 
trend. In NGC 3489, the ratio increases from the center out- 
wards. The [NIIl/Ho- ratio decreases very steeply in IC 1459, 
IC 4296, NGC 3607, NGC 4374, NGC 4697, NGC 5090, and 
NGC 5266. We show in Fig. [8] the individual BPT diagrams in 
different annuli for some galaxies of the sample. The [NIII/Hq- 
and [OIII]/HyS ratios in the four annuli for the whole sample are 
shown in Appendix B. 

Summarizing, the hardness of the ionizing continuum de- 
creases with galacto-centric distance. 

4.3. Central correlations with galaxy properties 

We investigated possible correlations between the emission line 
properties and the galaxy properties. Our results are summarized 
in Fig. m 

In the top panels of Fig.|9]we plotted the central (r < r^/lG) 
Ha and [NII]/16584 equivalent widths versus the central velocity 
dispersion cr^., the age, the metallicity, and the a/Fe enhance- 
ment. The Ha and [Nil] EWs quantify the strength of the emis- 
sion with respect to the stellar continuum. However, the [Nil] 
line is much less affected than the Ha by uncertainties or sys- 
tematics due to the stellar continuum subtraction. The scatter in 
the Ha and and [Nil] EW increases with cr^.. This still holds 
if the [OIII] emission is considered. If we exclude the Seyfert 



galaxy IC 5063, strong emission is observed only for galaxies 
with > 200 km s"' . The galaxies with the strongest emissions 
tend to be fast rotators. Also, all the low cr^ galaxies are fast ro- 
tators. Concerning the stellar population parameters, the trend 
is preserved only for the [Nil] line versus the [a/Fe] ratio (it is 
barely present in the Ha vs [a/Fe] plane). We notice that this 
trend could be driven by the contamination from the [NI]5199 
line to the Mgb index, implying larger derived [a/Fe] values for 
galaxies with stronger emission (Goudfrooij & Emsellem 1996Jl. 
The lack of a trend between the EW and metallicity is a bit sur- 
prising at the light of the well know metallicity-cr relation for 
ETGs. 

In the third row panels of Fig. |9] we show the central (r < 
re/16) [NII]/15684/Ha ratio, which is a measure of the hardness 
of the ionizing spectrum. In examining possible correlations be- 
tween this emission line ratio and the galaxy properties, we sep- 
arated the sample in high emission galaxies (H) and weak emis- 
sion ones (W) (see Section 3.3 and classification in Table 4). In 
fact, the determination of the Ha emission can be strongly af- 
fected by systematics in the underlying continuum subtraction 
in the W galaxies. To investigate the presence of correlations it 
is thus safer to consider only the H subsample. The first panel in 
the third row shows that there is a positive correlation between 
the central [NII]/Ha and cr^. The Spearman correlation coeffi- 
cient r, is 0.31 (with 36 data points), implying that there is a 6% 
probability that the variables are not correlated for a two-tailed 
test. However, if we restrict the analysis to the galaxies with 
EW([NII])> 3, the probability that there is no correlation drops 
to 4%. This result is in agreement with Phillips et al.( 1 19861 1, 
who derived a positive trend between the [Nil]/ Ha ratio and the 
galaxy absolute magnitude. On the other hand. Spearman tests 
show that there is no strong correlation between [Nil]/ Ha and 
the age, the metallicity, or the a/Fe enhancement of the galaxy 
stellar population. 

In the forth row panels of Fig. |9] we show the central 
[OIII]/H/? ratio, which is related to the ionization parameter. No 
clear trend is observed with 0"^, age, metallicity, and [a/Fe]. We 
only notice that the largest [OIII]/H/3 values are associated to the 
youngest galaxy ages. 
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Fig. 12. Oxygen abundance at r <ie/l6 obtained with the second 
calibration of Storchi-Bergmann et al. (JJ)98 ) against the values 
obtained with the first calibration. For comparison, the one-to- 
one relation is drawn as a solid line, while the dashed line indi- 
cates a systematic shift of 0.09 dex. 



The relation between the nuclear emission spectral class 
(Seyfert, LINERs, composite), and the stellar population prop- 
erties is show in Figs. [TO] and [TT] Seyfert galaxies tend to have 
ages younger than ~ 4 Gyr "Composite" galaxies have ages 
older than 5 Gyr. LINERs cover the whole age range from a few 
Gyrs to a Hubble time. For the three classes, the average ages 
are 3.8, 7.5, and 9 Gyr, respectively. The average metallicities 
are log Z/Zq ~ 0.39, 0.29, and 0.22, respectively. No significant 
difference is observed in the [ff/Fe] ratio. 

It emerges that Seyferts tend to have luminosity-weighted 
ages significantly younger than LINERs, likely because of recent 
star formation episodes. The only caveat in this statement is that 
the only bona fide Seyfert in our sample (IC 5063), does not have 
an age determination. 



5. Oxygen abundance 

If the ionizing source is known, one can derive the element abun- 
dance of the warm gas. In their nuclear (r < re/16) region, the 
majority of the galaxies in our sample are classified as LINERs, 
a few are classified as "Composites", and only IC 5063 is clas- 
sified as a bona fide Seyfert galaxy. We have also noticed a pro- 
gressive increase of the fraction of "Composites" with respect to 
pure LINERs moving from the galaxy center to more external 
regions. As discussed in Section 1, the ionizing mechanism in 
LINERs/Composites is still far from being understood. For this 
reason, we derived the oxygen metallicity for the R05-I-A06 sam- 
ple in the two extreme assumptions of 1) photoionization by hot 
stars, and 2) AGN excitation. A calibration in the case of shock 
heating is not present in the literature. 




[0/H] Hll 

Fig. 13. Top panel: oxygen abundance distri- 
butions ([O/H] = log(0/H) - log(0/H)o^ where 
12 + log(O/H)0 = 8.76 from Caff'au et al. l2008l) derived 
for the R05-HA06 sample at r <r(,/16. The dashed histogram 
was obtained with the Kobulnicky et al. ( 119991 ) calibration valid 
for HII regions, while the full histogram was obtained with the 
Storchi-Bergmann et al. (11998b calibration for AGNs. Vertical 
arrows indicate the average [O/H] values obtained with the 
two methods. Bottom panel: comparison of the metallicities 
obtained with the two methods for 20 galaxies in common. The 
solid line is the one-to-one relation. 



5.1. HII region calibration 

In HII region, the oxygen metallicity is determined through 
the standard R23 = ([OII]^3727 + [OIII]/l/}4959, 5007)/H;8 
parametrization, first introduced by Pagel et al. (1979), and 
then revised by other authors (Edmund & Pagel il984i McCall, 
Rybski & Shields 1985; Dopita & Evans [T986l McGaugh [T99Tl 
Zaritsky et al. 1994). The main problem with this calibration is 
that the relationship between R23 and the oxygen abundance is 
double valued, requiring some a priori knowledge of a galaxy's 
metallicity in order to determine its correct location in the 
upper or lower branch of the curve. A good discriminator is 
the [OIII]^5007/ [NII]/16584 ratio, which is usually less then 
K 100 for galaxies with 12 + log(0 /H) > 8.3, on the metal 
rich branch (Edmunds & Pagel Kobulnicky et al. [T999] l. 

The various calibrations, based on different photoionization 
models and HII region data, show a dispersion of ~ 0.2 dex 
in the metal rich branch. Kobulnicky et al. ( 119991 1 reports 
analytic expressions based on a set of photoionization models 
that fit the empirical calibrations within less than 0.05 dex, 
and that include the effect of the ionization parameter through 
O32 = [OIII]^^4959,5007/[OII]^3727. 

We computed the oxygen abundance for our sample adopt- 
ing the calibration in Kobulnicky et al. (119991 ) (hereafter Kob99). 
Because of the lack of star forming regions in ETGs, this cali- 
bration is valid if: 1) the gas is photoionized by hot old stars, 
likely PAGB stars, and 2) the calibration is still valid even though 
PAGB stars have different spectral shapes than high-mass main 
sequence stars. The same approach was adopted by Athey & 
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Fig. 14. Oxygen abundance distributions in the three more in- 
ternal annuli (r < r(,/16, rf/16 < r <xj%, re/8< r <re/4) obtained 
with the Kob99 calibration (18 galaxies, top panel), and with the 
SB98 calibration (6 galaxies, bottom panel). Vertical aiTows in- 
dicate the average values for the distributions. Notice that the 
metallicity increases from the galaxy center outwards with the 
Kob99 calibration, while the opposite trend is obtained with the 
SB98 caUbration. 



Bregman (I2009I I to determine oxygen abundances in a sample 
of 7 ETGs. We computed the R23 and O32 values using the emis- 
sion line fluxes in Table 2, and correcting for reddening through 
the values in Col 3 of Table 2. The derived oxygen metallicities 
at r <rp/16 are given in Col. 2 of Table 6. The errors were com- 
puted by combining the intrinsic uncertainty in the calibration 
(~ 0.15 dex, see Kob99), with the uncertainty in the emission 
line fluxes. We could not determine the oxygen abundance in 13 
galaxies, namely IC 5063, NGC 777, NGC 1297, NGC 1521, 
NGC 1947, NGC 2749 , NGC 3136, NGC 3489, NGC 3557, 
NGC 5011, NGC 6776, NGC 6958, and NGC 7007, since R23 
is larger than 10, and falls outside the valid calibration range for 
HII regions. 

There are 6 galaxies in common with the sample of Athey 
and Bregman (2009): NGC 3489, NGC 3607, NGC 4374, 
NGC 4636, NGC 5044, and NGC 5846. Excluding NGC 3489, 
we derived for the other galaxies central log Rx>, values of 0.50 
+ 0.15, 0.58 ± 0.11, 0.63 ± 0.1, 0.70 ± 0.05, 0.51 + 0.11. 



3 0.09 



0-58T,0'0.68V„, 



For the same galaxies, they derive 0.93 q j 
0.75"fog, 0.68"-^^o9. With the exception of NGC 3607, our and 
their results are in agreement within the errors. 



The oxygen metallicity distribution is shown in the top panel 
of Fig. [13] . The [O/H] = log(0/H) - log(0/H)Q values were 
computed adopting a solar abundance of 12 + log(0/H) = 8.76 
from CafFau et al. (2008). The average metallicity is 
[O/H] - -0.06. Some values are as low as ■x, 0.25 times solar, 
while the largest derived metallicity is « twice solar. 




[O/H] stars 

Fig. 15. Nebular oxygen abundance versus star oxygen abun- 
dance. The latter was derived from the Z and [a/Fe] values ob- 
tained in Paper III. In the top and bottom panels, the nebular 
metallicities were obtained with the Kob99 and SB98 calibra- 
tion, respectively. 



5.2. AGN calibration 



Storchi-Bergmann et al. d 19981 1 (hereafter SB98) derived oxy- 
gen metallicity calibrations from a grid of models assuming 
photoionization by a typical AGN continuum (the segmented 
power law of Mathews & Ferland 1987). The first calibra- 
tion (hereafter SBl) is in terms of [OIII]/1/14959,5007/ H/? 
and [NII]^^6548,6584/ Ha (Eq. (2) in SB98), while the sec- 
ond calibration (hereafter SB2) is in terms of [OII]/13727/ 
[OIII] ^4959,5007 and [NII]^^6548,6584/Hq' (Eq. (3) in SB98). 
The fitted calibrations are within ~ 0.05 dex of the models. 
The models were computed assuming a gas density ng = 300 
cm"^. The correction due to deviations from this density is 
-0.1 X log(ng/300), valid for 100 cm^^ < n^ < 10,000 cm^l 

We derived oxygen metallicities for the R05-I-A06 sample 
using both SBl and SB2. The electron densities n^, were calcu- 
lated from the [SII]/1/16717, 31 line ratios using the temden task 
within IRAF/STSDAS (Shaw & Dufour 1995), and assuming an 
electronic temperature Tg =10,000 K. This was possible only for 
some galaxies of our sample (see Table 5), because the quality 
of our spectra gets significantly worse beyond ~ 6600 A, and the 
presence of telluric absorption lines around 6800 A prevents in 
many cases a reliable fit to the [SII] doublet. The derived electron 
densities are given in Table 5. As evident from the [SII] fluxes 
given in Table 2, both increasing and decreasing density profiles 
are present within our sample. The oxygen metallicities were 
then corrected for the density dependence assuming Ug ^ n^. 
Because of the limited range of validity for the density correc- 
tion, we discarded the galaxies with n,. < 100 cm"^. The [O/H] 
values obtained with SBl and SB2 are compared in Fig. [12] The 
metallicities obtained with SB2 are on average ~ 0.09 dex larger 
than those obtained with SBl, in agreement with SB98 who de- 
rived a shift of ~ 0.11 dex between the two calibrations. The 
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Table 5. Electron densities n^, derived from [SII]6716/6731 as- 
suming T= 1 0,000 K 



Galaxy [Sll]6717/6731 n, (cm^ 



NGC 1052 


1.11 


+ 


0.03 


384.1 


NGC 1209 


1.19 


+ 


0.13 


254.8 


NGC 1297 


1.08 


+ 


0.10 


423.6 


NGC 1380 


0.98 


+ 


0.15 


653.9 


NGC 1453 


1.28 


+ 


0.09 


141.5 


NGC 1533 


1.22 


+ 


0.16 


212.0 


NGC 1553 


1.21 


+ 


0.30 


229.3 


NGC 1947 


1.32 


+ 


0.04 


95.0 


NGC 2749 


1.18 


+ 


0.08 


267.7 


NGC 2974 


1.10 


+ 


0.05 


385.1 


NGC 3136 


1.60 


+ 


0.11 


<26 


NGC 3268 


1.32 


+ 


0.13 


103.4 


NGC 3489 


1.21 


+ 


0.14 


223.1 


NGC 3607 


1.24 


+ 


0.20 


185.4 


NGC 3962 


1.27 


+ 


0.06 


147.1 


NGC 4374 


1.45 


+ 


0.17 


<26 


NGC 4636 


1.45 


+ 


0.15 


<26 


NGC 4696 


1.31 


+ 


0.08 


110.0 


NGC 5044 


1.51 


+ 


0.06 


<26 


NGC 5077 


1.26 


+ 


0.06 


158.7 


NGC 5090 


1.48 


+ 


0.20 


<26 


NGC 5193 


1.57 


+ 


0.50 


<26 


NGC 5266 


1.06 


+ 


0.04 


476.7 


NGC 5363 


1.27 


+ 


0.05 


148.8 


NGC 5813 


1.34 


+ 


0.22 


83.4 


NGC 5846 


1.10 


± 


0.21 


395.6 


NGC 5898 


1.55 


+ 


0.34 


<26 


NGC 6758 


1.42 


+ 


0.17 


<26 


NGC 6868 


1.43 


+ 


0.08 


<26 


NGC 6958 


0.77 


+ 


0.03 


1605.9 


NGC 7007 


1.53 


+ 


0.15 


<26 


NGC 7079 


0.95 


+ 


0.21 


732.3 


NGC 7097 


1.23 


+ 


0.06 


200.4 


NGC 7377 


1.43 


+ 


0.14 


<26 


IC 1459 


0.97 


+ 


0.04 


674.2 


IC 3370 


1.48 


+ 


0.05 


<26 


IC 4296 


1.10 


± 


0.15 


400.7 


IC5063 


1.08 


+ 


0.01 


433.3 



final metallicities are computed as the average between the two 
values. The results are given in Col. 3 of Table 6. 

We compare the results obtained with the Kob99 and with the 
SB98 calibrations in Fig. [13] The SB98 calibration provides an 
average oxygen abundance ~ 0.04 dex higher than that obtained 
with Kob99. Also, Kob99 provides a broader metallicity range 
than SB98. 

5.3. Abundance gradients 

We studied the behavior of the oxygen metallicity as function of 
galacto-centric distance. The [O/H] distributions obtained with 
Kob99 and SB98 in the three more internal annuli (r < rp/16, 
ie/l6 < I <r(,/8, re/8< r <tJ4) are shown in Fig. [141 The samples 
are not very large because of the increased difficulty in measur- 
ing the emission lines in the more external annuli. Emission lines 
are too faint to allow a reliable abundance measurement in the 
more external annulus (r(,/4< r <rg/2). The results obtained with 
the Kob99 and SB98 calibrations are contradictory: in the first 
case, the oxygen metallicity tends to increase from the galaxy 
center outwards. In the second case, the opposite behavior is 
observed. The latter trend is more consistent with the negative 
metallicity gradients of the stellar populations. 



Table 6. [O/H] values obtained the with Kobulnicky et al. (119991) 
calibration (Kob99), and with the Storchi-Bergmann et al. (119981) 
calibration (SB98). 



Galaxy 


log(0/H) + 12 (Kob99) 


log(0/H) + 12 (SB98) 


NGC 128 


8.44 ± 0.24 


— 


NGC 1052 


8.52 ±0.15 


8.67 ± 0.03 


NGC 1209 


8.63 ±0.19 


8.78 ± 0.04 


NGC 1297 


— 


8.59 ± 0.04 


NGC 1380 


8.81 ±0.20 


8.75 ± 0.04 


NGC 1453 


8.55 ±0.19 


8.76 ± 0.04 


NGC 1533 


8.76 ±0.18 


8.74 ± 0.04 


NGC 1553 


8.48 ±0.31 


8.87 ± 0.04 


NGC 2749 


— 


8.7 ± 0.04 


NGC 2911 


8.56 ±0.16 


- 


NGC 2974 


8.53 ±0.18 


8.86 ± 0.04 


NGC 3258 


8.94 ±0.17 


_ 


NGC 3268 


8.83 ±0.16 


8.85 ± 0.04 


NGC 3607 


8.86 ±0.18 


8.94 ± 0.04 


NGC 3962 


8.62 ±0.18 


8.83 ± 0.04 


NGC 4374 


8.79 ±0.17 


— 


NGC 4552 


8.98 ±0.16 


_ 


NGC 4636 


8.75 ±0.17 


_ 


NGC 4696 


8.71 ±0.16 


8.82 ± 0.04 


NGC 4697 


8.85 ± 0.26 


- 


NGC 5044 


8.66 ±0.16 


— 


NGC 5077 


8.40 ±0.17 


8.69 ± 0.04 


NGC 5090 


8.93 ±0.16 


_ 


NGC 5193 


9.01 ±0.16 


_ 


NGC 5266 


8.20 ± 0.22 


8.72 ± 0.04 


NGC 5328 


8.96 ± 0.23 


8.54 ±0.17 


NGC 5363 


8.32 ±0.21 


8.68 ± 0.03 


NGC 5813 


8.84 ±0.17 


— 


NGC 5846 


8.84 ±0.17 


8.71 ± 0.04 


NGC 5898 


8.75 ± 0.21 


— 


NGC 6721 


8.88 ±0.16 


_ 


NGC 6758 


8.63 ±0.18 


— 


NGC 6776 


_ 


_ 


NGC 6868 


8.48 ±0.19 


— 


NGC 6876 


8.97 ±0.16 


— 


NGC 6958 


- 


8.82 ± 0.04 


NGC 7079 


8.72 ± 0.28 


8.78 ± 0.05 


NGC 7097 


8.59 ±0.17 


8.67 ± 0.03 


NGC 7135 


8.64 ±0.16 




NGC 7192 


8.83 ±0.19 




NGC 7377 


8.60 ± 0.25 




IC 1459 


8.54 ±0.16 


8.78 ± 0.03 


IC 2006 


8.85 ± 0.21 




IC 3370 


8.49 ± 0.21 




IC 4296 


8.85 ±0.17 


8.88 ± 0.04 


IC5063 




8.69 ± 0.03 



5.4. Comparison witli stellar metallicities 

To understand the origin of the warm gas in ETGs it is important 
to compare the nebular metallicities with those of the stellar pop- 
ulations. We akeady derived the (total) stellar metallicities (Z) 
and the [a/Fe] ratios for the R05-I-A06 sample through the Lick 
indices in Paper III. The comparison with the nebular metallic- 
ities is not straightforward, since the emission lines provide a 
direct measure of the oxygen abundance, while the Lick indices 
depend on a mixture of different elements. In this paper, we com- 
pute the oxygen stellar metallicities from the Z and [a/Fe] values 
derived in Paper III, assuming that O belongs to the a-enhanced 
group (Ne, Na, Mg, Si, S, Ca, Ti, and also N). The compari- 
son between gas and star metallicities is shown in Fig. [15] The 
clear result is that, irrespective of the adopted calibration (Kob99 
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Fig. 16. Oxygen abundances derived with the Kob99 (top) and 
SB98 (bottom) calibration at r <rf/16 versus the galaxy central 
stellar velocity dispersion. Full symbols are the high emission 
galaxies (H), as defined in Table 4. The solid lines are the least- 
square fits to the data. The dashed and dotted lines are the re- 
lations derived by Thomas et al. ( 12005b and AO? for the stars, 
rescaled to account for the solar metallicity Zq = 00156 adopted 
in this paper. 



or SB98), the gas metallicity tends to be lower than the stellar 
metallicity. The effect is more severe for the galaxies with the 
largest stellar metallicities. A Spearman test shows that in both 
cases the star and gas metallicities are poorly correlated. 

We plot the central (r <ijl6) nebular metallicities obtained 
with Kob99 and SB98 versus the galaxy stellar velocity disper- 
sions in Fig. [16] Linear square fits provide an increasing trend of 
[O/H] with (Tc in both cases. However, the Spearman correlation 
coefficients are low (r, = 0.25 with N = 39 degrees of freedom 
for Kob99, top panel, and r., = 0.08 with N - 22 degrees of 
freedom for SB98, bottom panel), indicating that the correla- 
tions are weak. For comparison, we plotted in Fig.[T6]the stellar 
metallicity-cTc relations derived by us in Paper III (A07) and by 
Thomas et al. ( 120051 ) for low density environment ETGs. The 
A07 relation was appropriately re-scaled to account for the solar 
metallicity adopted in this paper Q In Fig. [161 the fits provide 
[0/H]g,s ~ -0.08 and [0/H]gas ~ -0.03 at cr, = 200 km s ' for 
Kob99 and SB98, respectively. For comparison, the A07 relation 
gives a stellar metallicity of ~ 0.18 dex at the same cr, i.e. 0.26 
dex and 0.21 dex larger than those obtained for the warm gas 
with Kob99 and SB98. The discrepancy is even larger if we con- 
sider the Thomas et al. (120051) relation, which provides a metal- 
licity of 0.34 dex at cr^ = 200 km s"', i.e. ~ 0.4 dex higher than 
the nebular metallicity. We notice that the discrepancy between 
gas and stars persists (and is even more significant) if we con- 
sider only the high emission (H) galaxies (see classification in 
Table 4). 



' Zo = 0.018 in A07, while we adopt Zq = 0.0156 from 
Caffau et al. (2009 1 in this paper. This implies a correction term of 
log(0.0 18/0.0 156) = 0.06 dex to be applied to the A07 relation. 



♦ • X/radio CS" 

2 Zs # EW(Ha) <1 - 




-1 -0.5 0.5 



log([NIl]/Ha) 

Fig. 17. BPT diagram for our galaxies at r <re/16 with superim- 
posed the shock models of Allen et al. I2008l and the dusty- AGN 
models of Groves et al. (2004). Full dots indicate objects where 
compact nuclear X-ray and/or radio sources have been detected 
(X/radio CS), consistently with the presence of an AGN, while 
stars are for objects with equivalent width in Wa < 1 A , com- 
patible with ionization from PAGBs. Notice that three galaxies 
(NGC 777, NGC 4552, and NGC 3557) are X/radio CS and at 
the same time have EW(Hq'< 1). The shock models have solar 
metallicity, densities from n - 0.01 cm"-' to 1000 cm"-', veloci- 
ties from V = 1 00 km ' to 500 km 5" ' , and magnetic field B - \ 
jiG. The plotted AGN models have solar and twice solar metal- 
licities, n - 1000 cm"^, and ionization parameters from logU ~ 
-4 to 0. The models were downloaded from the Brent Groves 
web page http://www.strw.leidenuniv.nl/~brent/itera.html using 
the ITER A tool. " ~~ 

6. Comparison with IVIodels 

The ionizing source in LINERs is still not well understood. 
Among the proposed mechanisms, we recall sub-Eddington ac- 
cretion into a super massive black hole (e.g., Ho |1999bl Kewley 
et al. 120061 Ho |2009a| l, photoionization by old PAG B stars 
(e.g. Trinchieri & di Serego Alighieri 1991 ; Binette et al. 19941 
Stasiriska et al. .2008 J , fast shocks ( Koski & Osterbrock , 19761 
Heckman [19801 Dopita & Sutherland [T995l Allen et al. l2008l l. 

The first hypothesis is strongly supported by our current 
knowledge of the demography of central BHs based on di- 
rect dynamical searches (Magorrian et al. 119981 Ho II999at 
KormendV i2004 ). according to which massive BHs appear to be a 
generic component of galaxies with a bulge. As shown by Ho et 
al. ( 1997c), local stellar mass loss can supply the fuel necessary 
to accretion rates of M = 0.001 - QAM^yr^^ . Thus, radiatively 
inefficient accretion flows onto a central BH (see Ouataert |200II 
for a review) provide an attractive solution to explain the ori- 
gin of LINERs. Additional evidences for the presence of AGN 
activity come from X-ray and radio observations, which have re- 
vealed the presence of compact nuclear X-ray and radio sources, 
sometimes with detected variability, in a significant fraction of 
LINERs (see e.g. Terashima ,1999i Flohic et al. ,2006i Filho et 
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Fig. 18. Distribution of L[oin]/cr'* ratio for Seyferts, Composites and LINERs (left panel), and for the total sample and galaxies 
classified as AGN and non AGN from X-ray/radio data (right panel). 



al. 120061 Gonzalez-Martin et al. l2009l Plan et al. lIOTOl l. On 
the other hand, some LINERs lack evidence of a central AGN, 
leaving space to alternative explanations for the observed emis- 
sion lines. X-ray and radio studies are available in the litera- 
ture for a few galaxies in our sample. Compact nuclear X-ray 
and/or radio sources have been detected in IC 5063 (Koyama 
et al. [T992] |. NGC 777 (Ho & Ulvestad [2 00D, NGC 1052, 
NGC 4374, IC 1459 (Gonzalez-Martin et al. '2009^, NGC 5363 
(Gonzalez-Mart in et al. l2009l Filho et al. 2006), NGC 1553 
(Blanton et al. [200Tl l. NGC 3 557 (B almaverde et al. 2006), 
NGC 5090 (Griitzbauch et al. 120071), IC 4296 (Pellegrini et 
al.|2003 ), NGC 5077, (Filho et al. l2006] l. and NGC 4552 (Nagar 
et al. i2002i Filho et al. |2004] |. Other galaxies observed in X- 
ray and/or radio lack evidence for AGN activity; NGC 3607, 
N GC 46 36, NGC 5813, (Gonzalez-Martin et al. l2009l Filho e t 
al. 12006]) . NGC 4696, NGC 5 846 (G onzalez-Martin et al. l2009] l. 
NGC 6876 (Machacek et al. l2003T l. High quality X-ray imag- 
ing has revealed in NGC 4636 the signature of shocks, probably 
driven by energy deposited off'-center by jets (Jones et al. i2002i 
Baldi et al. l20^ . 

Shock heating was initially proposed as a viable excita- 
tion mechanism to solve the so-called "temperature problem" of 
Seyferts and LINERs, in which the electron temperatures ob- 
served were found to be systematically higher than predicted 
by photoionization models. Indeed, Koski & Osterbrock d 19761 ) 
and Fosbury et al. (il978l ) argued that shocks are essential to ex- 
plain the observed [OIII] line ratios in the prototypical LINER 
NGC 1052. Dopita & Sutherland ( [1995] ) and Allen et al. (12008] ) 
presented extensive grids of high velocity shock models, and 
showed that the emission line ratios observed in LINERs can be 
modeled in terms of fast shocks ( 1 50 - 500 km s" ' ) in a relatively 
gas poor environment. 

The findings that the emission-line flux correlates with the 
host galaxy stellar luminosity within the emission-line region 
(Macchetto et al.ll996), and that the line flux distribution closely 
follows that of the stellar continuum (Sarzi et al. [2006 2009), 
have been presented as evidences in support of the PAGE sce- 
nario. These results suggest in fact that the sources of ionizing 
photons are distributed in the same way as the stellar popula- 
tion. Recendy, Sarzi et al. (I2009i ) tried to model the H/? EW 



spatial distribution under the assumption that the ionizing pho- 
tons originate from an underlying old stellar population. They 
showed that to obtain an almost constant EW it is necessary 
to assume that the gas radial profile decrease more gently than 
the ionizing radiation, or alternatively that the ionizing photons 
are more concentrated than the bulk of the stellar populations. 
Indeed, the presence of negative metalhcity gradients in ETGs 
(e. g., Annibali et al. 120071 ) could explain this eff'ect, because in 
general the ionizing flux from old populations is larger for higher 
metallicities. The Stasihska et al. (12008b models are able to re- 
produce the emission line ratios of LINERs provided that the 
ionization parameter is sufficiently high (logU > -4). Binette 
et al. (1994) and Cid Fernandes et al. ( 120091 ) demonstrated that 
photoionization by PAGE stars can only explain LINERs with 
relatively weak emission lines. Einette et al. (11994 ) estimated 
that the Ha equivalent width produced by a 8 Gyr old and a 13 
Gyr old simple stellar populations amounts to 0.6 A and 1.7 A, 
respectively. Stronger emission lines require far more ionizing 
photons than old stars can provide. 

A first-order test on the validity of the PAGE scenario 
consists in a comparison of the observed and predicted EWs. 
Following the computations of Einette et al. ( 119941 ). we assume 
that photoionization by old stars produce a typical emission of 
1 A in Ha. Comparing this value with our measured central 
(r < rg/ 16) EWs, we obtain that PAGE photoionization alone 
can explain the observed nuclear emission in only 1 1 out of the 
49 galaxies classified as LINERs/Composites in our sample (i.e. 
in a; 22% of the LINERs/Comp sample). This implies that, for 
the majority of the cases, some mechanism other than PAGE 
stars must be at work in the central regions of LINERs. On the 
other hand, as we progressively move toward larger radii, an in- 
creasing number of LINERs enter the EW(Ha)~ 1 A transition 
region, where the emission strength is compatible with photoion- 
ization by PAGE stars alone. Thus, while our computations ex- 
clude that PAGE stars play a significant role in the central galaxy 
regions, we can not rule out that they are an important excitation 
source at larger radii. Interestingly, Sarzi et al. (2 009) quote a 3" 
radius for the maximum extent of central activity, which for their 
sample galaxies corresponds on average to 0. 1 1 r^,. 
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Fig. 19. Shock models of Allen et al (2008) compared with our 
r <re/16 Ha luminosities, corrected for reddening. The shock 
models are plotted for 4 density values (n - 0.01, 1, 10, 100 
cm~^), and magnetic fields B = 0. 1 (dotted line), B = 1 (solid 
line) and B = 10 fiG (dashed line). The velocity in the abscissa 
is the shock velocity for the models, and V3cr^- for the galaxies. 
The different symbols used to plot the data are the same as in 
Fig. [17] We indicate the name of the galaxies with the largest 
Lua values. NGC 777, once corrected for the high extinction 
(E(B-V)=1.6), has Lh^ ~ 0.06 ergs^'cm"^ and falls outside the 
plot boundaries. 

One of the arguments presented in support of the PAGE pho- 
toionization scenario is that the emission flux closely follows 
the stellar continuum. At odds with other results (e.g. Sarzi et 
al. I2009I I, we have shown in Section 3.3 that the EW (emission 
flux over continuum flux) decreases with radius. Is this incom- 
patible with the PAGE scenario? Answering this question is not 
straightforward. First of all, we have shown that the strongest 
decrease in EW occurs from the nuclear to the second (it,/l6 
< r <re/8) annulus, being much gentler at larger radii. Indeed, 
our data lack the depth necessary to reliably trace the emission 
strength beyond 1^/4-. Second, neither a well established decreas- 
ing trend in EW would be per se incompatible with the PAGE 
scenario: in fact, because of the presence of age/metallicity gra- 
dients in the stellar population, the ratio of the ionizing flux over 
the optical continuum is not necessarily constant. A reliable ra- 
dial prediction of the ionizing photon budget stems on a detailed 
analysis of the combined effect of age and metallicity on the 
PAGE population, which is beyond the aims of this paper 

Since the nuclear emission requires a mechanism differ- 
ent than PAGE photoionization, we compare in Fig. [T7]our (r 
<re/16) data, corrected for reddening, with the dusty AGN mod- 
els of Groves et al. (I2004I I. and with the fast shock models of 
Allen et al. (2008). Among the data points, we identify ob- 
jects with compact nuclear X-ray and/or radio source detections, 
likely hosting an AGN, and objects with equivalent width in 
Ha < 1 A , compatible with sole ionization from PAGEs. We no- 



tice that fliree galaxies (NGC 777, NGC 4552, and NGC 3557) 
present evidence for AGN activity from X ray or radio obser- 
vations, and at the same time have EW(Ha)< 1. Comparing 
dusty models with reddening-corrected data is not completely 
self-consistent; however the effect of extinction in the EPT dia- 
gram is negligible. 

Eoth the AGN and shock models were developed using 
the co de MA PPINGS III (Dopita et al. 1982, Sutherland & 
Dopita I1993I I, implemented with dust and radiation pressure 
(Groves et al. '2004). The AGN models were developed for dif- 
ferent metallicities, and assuming a power law with slope a be- 
tween -1.2 and -2.0 for the ionizing spectrum. For the AGN 
models in Fig.[T7] a = -1.4, Z - Zq and 2 Z©, uh =1000 cm"-', 
and log U between -4 and 0. LINERs and Seyferts are both 
explainable with being AGNs, provided that the formers have 
much lower ionization parameters than the latters. According to 
these models, the galaxies in our sample should have gas metal- 
licities much above solar This is in disagreement with the oxy- 
gen metallicities derived in Section 5 with the Storchi-Eergmann 
et al. ( 1998) calibration. The discrepancy arises from the use of 
a power law in the Groves et al. (2004) models compared to the 
typical AGN continuum of Mathews & Ferland (I1987I I adopted 
by SE98. SE98 noticed in fact that models obtained with power- 
law ionizing continua yield systematically larger (up to 0.5 dex) 
metallicities than models that adopt a typical AGN continuum. 

Under the assumption that all the central LINERs activity 
is powered by a sub-Eddington accretion AGN, we can attempt 
a determination of the accretion rate from the L[oiii]/o-'* ratio 
(where L[oiii] is the extinction corrected [OIII] luminosity), pro- 
portional to L/Ledd (Heckman et al. 2004). The result is shown 
in Fig. [18] The accretion rate increases systematically along 
the sequence Composite — > LINERs Seyferts (left panel) in 
agreement with recent findings by Ho (2009a1). Also, galaxies 
with compact nuclear X-ray/ radio sources, likely hosting an 
AGN, tend to have larger accretion rates than non -AGN galaxies 
(right panel). Notice that since L oc tjLedd ^ tjct^, and tj « 1 
in the sub-Eddington accretion regime, the AGN luminosity can 
vary significantly for a fixed cr depending on the accretion rj. 
This implies that it is necessary to assume a maximum, sub- 
Eddington, accretion rate to reproduce the observed trend be- 
tween the emission strength and cr. 

The Allen et al. (I2008I I shock models are dust free, and were 
computed for different chemical abundances, pre-shock densi- 
ties n from 0.01 to 1000 cm"-', velocities up to 1000 km s"', and 
magnetic parameters Bln^l^ from lO""* to 10 jjG cm-'^-, where E 
is the pre-shock transverse magnetic field. In Fig. [T7] we plot- 
ted the models with solar metallicity and E= 1. With these pa- 
rameters, the bulk of LINERs in our sample is well reproduced 
by pre-shock densities n - 0.01 - 10 cm"^, and shock veloc- 
ities V = 200 - 500 km i"'. Using a higher magnetic field 
implies larger densities for our sample. For instance, for E=10 
I^G, our data are consistent with pre-shock densities as high as 
n X 100 cm"^. This range in E is typical of the ISM (e.g.. Rand 
& Kulkai-ni 1989 ). 

Are these densities consistent with those derived from the 
observed [SII](6717)/[SII](6731) ratios? As shown in Table 5, 
the majority of our galaxies have ratios > 1, implying ng < 
600 cm"^ for an assumed T= 10,000 K. A significant fraction of 
galaxies have ratios around 1.45, which is the saturation limit 
for n— > 0. These values are consistent with pre-shock densi- 
ties as low as n<10 cm""*. Indeed, as the density of the ma- 
terial passing through the shock increases, the density indi- 
cated by the [SII] ratio also increases. However, the amount 
of compression depends critically upon the magnetic parameter 
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(Dopita & Sutherland [T995T l. In the Allen et al. (l2008l l models, 
[SII](6717)/[SII](6731) ~ 0.7 and ~ 1 for a pre-shock density 
n=10 cm"^, and B=l and 10 fiG, respectively. 

In the fast shock scenario, where does the mechanical en- 
ergy come from? Dopita & Sutherland (1995) distinguish be- 
tween two possible types of gas dynamical flows: jet-driven and 
accretion-driven flows. Evidence for jets in LINERs comes from 
observations revealing the presence of low-luminosity AGNs as- 
sociated with radio sources with a high brightness temperature 
and/or elongated radio morphology (e.g., Nagar et al. 1200 11 1. In 
accretion-driven flows, the shock is either the result of the col- 
lision of cool clouds which have dropped out of a cooling flow 
into the potential of the galaxy (e.g. Crawford & Fabian [1992]l 
or have been accreted from a companion galaxy. Shocks may 
also occur in the outer accretion disk about the central object 
either trough the dissipation of turbulence, dynamical intsabili- 
ties, or spiral shocks, (e.g. Dopita et al. |1997l l. We suppose that 
the mechanical energy could also come from turbulent motion of 
gas clouds, injected by stellar mass loss or accreted, within the 
potential well of the galaxy. Fig.fTTlshows that the shock veloci- 
ties required to reproduce the distribution of LINERs in the BPT 
diagram range between 200 and 500 km/s. Are the mentioned 
mechanisms able to account for such high shock velocities? It 
seems that jet driven flows or accretion into a central massive BH 
can weU do the job (e.g., Nagar et al. 2005 , Dopita et al. |1997! ). 
On the other hand, it is not straightforward to establish if, in 
absence of these sources, pure gravitational motion of clouds 
within the potential well of the galaxy can produce collision at 
the required velocities. Close to the center of a relaxed stellar 
system, the gas should quickly relax into a disk-like configura- 
tion and orbit at a velocity close to the circular one, even though 
not necessarily in the same direction of stars. Indeed, this is con- 
firmed by observations, showing regular gaseous disk structures 
and coherent gas motions (see e.g. Zeilinger et al. 1996, Sarzi 
et al. I20061 I. If most of the clouds rotate in the same direction, 
it would be then rather difficult for them to collide at the high 
speed required. On the other hand, observations detect an addi- 
tional gaseous component of high gas velocity dispersion, in the 
range 150-250 km/s (Bertola et al [19951 Zeilinger et al. [19961 
Emsellem et al 2003). Sarzi et al. (I2006I I found that cr^„5 is gen- 
erally smaller than the stellar velocity dispersion cr», but in some 
cases (Tgas ~ cr», either only in the central regions or over most 
of the field. For a sample of 345 galaxies, with Hubble type from 
E to Sbc, Ho 2009b found that the gas dispersions strongly cor- 
relate with the stellar dispersions over the velocity range of cr a; 
30-350 km s"', such that o-gas /cr« x 0.6 - 1.4, with an average 
value of 0.80. All these observations suggest that random mo- 
tions are crucial for the dynamical support of the gas. Projected 
velocity dispersions in the range 150-250 km/s translate, under 
the assumption of isotropy, into intrinsic 3D velocity dispersions 
in the range 260-430 km/s. This is a crude estimate of the in- 
trinsic velocity dispersions since studies indicate deviation from 
isotropy (e.g. MeiTifield et al. 2001 ). Given these values, colli- 
sions at a speed between 200 and 500 km/s can not be excluded. 

We compare the shock models of Allen et al. ( I2008I I with our 
data in a luminosity vs velocity plane in Fig. [19] For the data, we 
adopted a velocity of V3cr^-, which is a first-order approximation 
for the shock velocity under the assumption that the shocks are 
due to motions of gas clouds in the galaxy potential well. The 
models are plotted for four density values (n = 0.01, 1, 10, 100 



^ The total 3D velocity dispersion satisfies cr- = a\+a^.+al, because 
of the Pythagoras theorem and basic statistics of the second moments. 
If the dispersion is isotropic, cr ., = cr,, = a- by definition, and cr^ = 3cr^. 



cm"^), and for three values of the transverse magnetic field 
{B - 0.1, 1, 10juG).The shock luminosity increases significantly 
with the shock velocity and with the gas density. The dependence 
on the magnetic parameter is almost null in Ha (and modest in 
[Nil]). The observed nuclear Ha luminosities were coiTected for 
reddening through the values derived in Table 2. Some galax- 
ies have very large reddening, associated to large uncertainties, 
and display quite high luminosities once corrected for extinc- 
tion (NGC 777, NGC 1521, NGC 7007, NGC 3557). Their Lh^ 
should be taken with caution. Fig. [T9]shows that, for almost all 
the galaxies in our sample, the observed luminosities are well 
reproduced by models with pre-shock densities n < 10 cm"^. 
This does not proof that shocks are the excitation mechanism in 
LINERs, but show that such a scenario is compatible with the 
data. We notice that shocks could provide a natural explanation 
for the observed half-cone-shaped distribution in the line lumi- 
nosity/EW vs cr plane (see also FigO. They could also explain 
the observed [NII]/Ha vs cTc coiTelation (FigO.The shock mod- 
els predict in fact an increase of the [NII]/Ha line ratio with 
the shock velocity. This is because higher velocity shocks result 
in harder and more luminous ionizing spectra. It is interesting 
that the weak emission galaxies in our sample, compatible with 
pure photoionization by PAGB stars at their center, are those 
deviating from the [NII]/Ha vs (Tc relation in Fig|9l Phillips et 
al. ( 1986) also found that the [NII]/Ha ratio increases with the 
galaxy absolute magnitude in a sample of LINERs. They inter- 
preted this behavior as due to differences in metallicities or to 
due to differing ionization conditions. 

7. Summary and Conclusions 

Optical emission lines are detected in 89% of our sample. This 
large fraction is somehow expected, since our sample was se- 
lected among a compilation of galaxies showing traces of ISM, 
and it is thus biased toward the presence of emission lines. This 
fraction is higher than that derived by Phillips et al. (11986!) (55%- 
60%), Macchetto et al. ([T996] l: (72% - 85%), Sar zi et al . ([20061) 
(75%), Yan et al. (TME) ( 52%), and Serra et al. ([2008l l (60%), 
for other samples of early-type galaxies. If we consider only the 
galaxies in our sample with EW(Hq' -h [NII]6584)>3 A, the de- 
tection fraction drops to 57%. 

The incidence and the strength of emission in our sample is 
independent on the elliptical or lenticular morphological classes. 
These results are in agreement with Phillips et al. ( 1986), while 
on the other hand Macchetto et al. ( 11996! ) and Sarzi et al. (2006| 
found a higher incidence of emission in lenticular than in elhpti- 
cal galaxies. Following the studies of Cappellari et al. ( !2007! l and 
Emsellem et al. (2007 ), we attempted a classification between 
fast and slow rotators in our sample. Fast and slow rotators ac- 
count for ~ 70% and ~ 30% of the sample, respectively. This is 
interesting given the fact that the sample is composed of ~ 70% 
ellipticals and ~ 30% lenticulars. Emission does not depend on 
the fast and slow rotator classification in our sample. 

The EW of the lines tends to decrease from the center 
outwards. The [Nil] line presents a steeper decrease than the 
Ha. However, given the uncertainties, we consider significant 
only the decrease from the center to ~ 0.1 r^.. Deeper observa- 
tions are needed to trace with higher confidence the EW out to 
larger radii. Some galaxies deviate from the general trend, ei- 
ther because they have a flatter emission profile (NGC 3489, 
NGC 7007, NGC 7377, NGC 5898), or because their emis- 
sion is nuclear concentrated (IC 4296, NGC 4374, NGC 5090). 
Previous narrowband imaging and integral field spectroscopy 
studies have revealed extended distribution for the ionized gas 
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(Caldwell fT984l Goudf rooij et al. ll994bt Macchetto et al. [19961 
Sarzi etal. '20061 120091) . 

Extinction was derived from the observed Hafli/3 flux ra- 
tio, assuming an intrinsic value of « 3.1 for AGN-like objects 
(Osterbrock 1989). This choice was driven by the fact that our 
galaxies show AGN-like [NII]/Hq' ratios. For the majority of 
the galaxies, the derived E(B - V) values are lower than 0.3. 
Some galaxies display however very large reddenings, up to 
E(B - V) ~ 1.5 or even more. Since the observed continuum is 
incompatible with such large values, our result suggests that in 
these galaxies the dust has a patchy distribution. This is consis- 
tent with narrowband imaging studies revealing the presence of 
dust lanes and patches in early-type galaxies (e.g. Goudfrooij et 
al. ll994bl) . 

We performed a spectral classification for our sample 
through the classical [OIII]/Hy6 versus [NII]/Ha' diagnostic di- 
agram (BPT diagram). The classification based on the central 
r<r(./16 emission lines is provided in Table 4. About 72% 
of the galaxies present a central LINERs emission activity. 
Seyferts account for ~ 9% of the emission sample. However, 
only IC 5063 can be classified as a bona fide Seyfert galaxy. 
The other 4 Seyferts (NGC 3489, NGC 777, NGC 7007, and 
NGC 6958), have very large errors in [OIII]/HyS (because of 
the faint emission in H/?), and the classification is less robust. 
We notice however that NGC 3489 and NGC 777 wer e classi- 
fied as Seyferts also by other authors (Ho et al. 11997b' Sarzi et 
al. |2006). 7 galaxies (~ 12% of the emission line sample) are 
"Composites" (or Transition objects): NGC 3258, NGC 4552, 
NGC 5193, NGC 5328, NGC 6721, NGC 6876, and IC 2006. 
They have [NII]/Hq' ratios intermediate between HII regions 
and LINERs/Seyferts. In our sample, these galaxies display the 
weakest emission lines. For the remaining 7% of the emission 
line galaxies we can not provide a classification, either because 
the lines are too weak, or because the spectra are too noisy. 

We derive spatial trends in the emission Une ratios: moving 
from the center toward annuli of increasing galacto-centric dis- 
tance, we detect a clear decrease in the [NIIl/Ho- ratio. The de- 
crease in [OIII]/H/3 is less clear It results that galaxies classified 
as LINERs from to their nuclear line ratios, shift toward the re- 
gion of "Composites" if one consider the emission from more 
external regions. This could be interpreted either with a decrease 
in the nebular metallicity, or with a progressive "softening" of 
the ionizing spectrum. 

We have investigated possible relations between the central 
line emission properties and the host galaxy properties. The main 
results are: 

1. Seyferts have young luminosity weighted ages (< 5 Gyr), 
and are on average younger than LINERs. LINERs span 
a wide age range, from a few Gyrs to a Hubble time. 
"Composite" have ages older than 5 Gyr, at variance with 
the idea that they may originate from a combined contribu- 
tion of AGN and star formation. 

2. Excluding the Seyferts, the spread in the emission equiva- 
lent width ([Nil], Ha or [OIII]) increases with the galaxy 
central velocity dispersion. Low-cr galaxies have weak emis- 
sion lines, while high-cr galaxies show both weak and strong 
emission lines. Equivalent widths in [NII]/16584 larger than 
~ 5 A are found only in galaxies with ctc > 200 km s"' . 

3. If we consider only the high emission subsample (H) (see 
Table 4), where the emission line ratios are more reliable, 
a positive correlation exists between the [NIIl/Ha ratio and 

(Tc. 



4. Excluding the Seyferts, which have larger [OIII]/HjS and 
younger ages, no relation between the emission line ra- 
tios and the galaxy age, metalhcity, or a/Fe enhancement is 
found. 

7. 1 . The AGN-starburst connection in Seyferts 

The first result is in agreement with previous findings that 
LINERs are older than Seyferts (e.g., Kewley et al. 20061), 
and supports the idea that the star formation and the AGN 
phenomena co-exist (Terlevich et al. 119901 Cid Fernandes & 
Terlevich 119951 Heckman et al. 119971 Gonzalez Delgado et 
al. 2001; Kauff^mann et al. TmJ: Cid Fernandes at al. [20051 
Davies et al. 2007; Riffel et al. 2009). The young luminosity 
weighted ages in the Seyferts of our sample are likely the result 
of recent star formation episodes superimposed to a several Gyrs 
old stellar population (see Paper III). Davies et al. (120071 ) provide 
evidence for a delay of 50- 1 00 Myr between the onset of star for- 
mation and the subsequent fueling of the black hole. Given the 
short duration of the AGN phenomenon (< 100 Myr), this sug- 
gests that star formation has occurred a few hundreds Myrs ago 
in the Seyfert galaxies of our sample. Because of the short dura- 
tion of the AGN phase, we also expect to observe galaxies with 
relatively young luminosity-weighted ages but no signatures of 
current Seyfert activity. This is in agreement with the fact that 
LINERs span a wide luminosity-weighted age range, from a few 
Gyrs up to a Hubble time. 

7.2. The ionization mechanism in LINERs 

Several studies have demonstrated that a minimum level of ion- 
izing photons are produced in any evolved stellar systems. Thus, 
if cold gas is present, some level of nebular emission is always 
expected in early-type galaxies. However, the critical question 
is: can PAGB stars provide all the ionization observed in early- 
type galaxies? In answering this question, we should first of all 
distinguish between the central LINERs activity, and the more 
extended LINERs-like emission. Indeed, through the BPT clas- 
sification carried out in annuli of increasing galacto-centric dis- 
tance, we have demonstrated that ETGs have both a central 
and an extended LINERs-like emission. When comparing the 
stronger nuclear r<r(./16 emission with the results of the pho- 
toionization models from Binette et al. (.1994 ), we obtain that 
only 1 1 out of the 49 galaxies classified as LINERs/Composites 
in our sample (i.e. in » 22% of the LINERs/Comp sample) can 
be explained with photoionization by old stars alone. Of these 
1 1 galaxies, 6 are classified as LINERs and 5 as Composites. 

For the other 78% of LINERs/Composites, some mechanism 
different than photoionization by PAGB stars must be at work. 
We warn however that this fraction is very uncertain, since there 
is no precise description yet for the evolution and number of 
PAGB stars (e.g.. Brown et al. [20081 ). 

On the other hand, we can not rule out the importance of hot 
old stars as a photoionization source at larger radii. Indeed, the 
emission strength progressively decreases toward larger galacto- 
centric distances, so that and an increasing number of LINERs 
become compatible with pure PAGB stars photoionization. 

Summarizing, we can not exclude a scenario in which more 
than a source of ionizing photons are present, with different roles 
at different radii. We can think of a transition region from the 
center, where PAGB stars fail in accounting alone for the ioniz- 
ing photon budget, to more extended regions, where they may do 



22 



Annibali et al.: Nearby early-type galaxies with ionized gas. 



well the job, in agreement with recent findings ( Sarzi et al. \2009\ 
Eracleous et al. \201Cf). 

Studies of LINERs based on large sample, such as the 
SDSS, have been recently presented (Stasinska et al. (120081 Cid 
Fernandes et al. (2009). However, we warn against direct com- 
parisons of our findings with the results from SDSS for at least 
two reasons: the 3"-wide aperture of the SDSS usually encom- 
passes large Kpc-scale regions where diffuse emission contribute 
to the nebular fluxes; and our sample is biased toward the pres- 
ence of emission lines. 

Once ascertained the limited role of PAGE stars in the nu- 
clear LINERs region, we investigated other possible scenarios, 
more specifically photoionization by sub-Eddington accretion 
rate AGN and fast shocks. Direct evidence for the presence of 
AGNs comes from the detection of compact nuclear X-ray/ ra- 
dio sources. High resolution X-ray and/or radio data are avail- 
able for 18 galaxies in our sample, 16 of which are classified as 
LINERs. Among the 16 LINERs, evidence of an AGN is found 
in 10. Even if the statistics are small, this result suggests the pres- 
ence of low accretion rate AGNs in ~ 62% of LINERs. However, 
the presence of an AGN does not exclude the coexistence of 
shock heating: indeed, jet driven flows and dissipative accre- 
tion disks in active galaxies with strongly sub-Eddington accre- 
tion have ben proposed as possible sources of shocks (Dopita & 
Sutherland (1995), Dopita et al. 119971 ). and bubbles likely pro- 
duced by an AGN outburst have been observed in LINERs (e.g. 
Baldi etal. 2009). 

Both AGN and shock models reproduce the central emission 
line ratios in our sample. In the dusty AGN models of Groves et 
al. (2004), LINERs require a ionization parameter logU < -3, 
significantly lower than for Seyferts. This is consistent with a 
scenario in which LINERs are AGNs accreting at a much lower 
rate than Seyferts. LINERs can also be modeled in terms of fast 
shocks (200-500 km s"') in a relatively gas poor environment 
(pre-shock densities n x 0.01 - 100 cm"^, depending on the mag- 
netic field), with gas metallicity ~ solar. These low pre-shock 
densities are consistent with the observed [SII]6717/[SII]6731 
ratios. A critical point concerns the high shock velocities re- 
quired to explain the emission line spectrum in LINERs. Jet 
driven flows or accretion into a massive BH can well do the 
job (Nagar et al. 120051 Dopita et al. II997I ). We suppose that 
an alternative source of mechanical energy comes from turbu- 
lent motion of gas clouds within the potential well of the galaxy. 
This scenario can not be excluded since several studies revealed 
a gaseous component of high velocity dispersion, in the range 
150- 250 Km/s, in the center of early-type galaxies (e.g. Bertola 
et al. [T995l Zeilinger et al.fT996', Sarzi et al. 12006] ). which trans- 
lates into higher intrinsic 3D velocity dispersions. This scenario 
is very appealing since it would naturally explain the observed 
trends between the emission strength, and the [NII]/Ha' ratio, 
with the galaxy stellar central velocity dispersion. 

A solid result from our study is that the [NII]/Hff ratio de- 
creases with increasing galacto-centric distance. This trend may 
reflect metallicity gradients in the gas, similarly to what ob- 
served for the stars. However, this is difficult to fit in a scenario 
in which the gas has an external origin (see Section 7.3). If not 
due to metallicity variations, the [NII]/Hq' gradient should re- 
flect variations in the properties of the ionization mechanism. 
In the shock scenario, this requires decreasing shock velocities 
with increasing galacto-centric distances. Unfortunately, we lack 
the spectral resolution necessary to determine the gas velocity 
dispersion and, ultimately, to ascertain the presence of a radial 
gradient in the shock velocity. 



A large fraction (40/65) of our sample was observed with 
Spitzer-IRS. In a forthcoming paper, we will use MIR ionic and 
molecular emission lines, and Polycyclic Aromatic Hydrocarbon 
features, to further investigate the powering mechanism in 
LINERs, and better disentangle the contributions from AGN 
photoionization and shock heating. 

7.3. Origin of tiie gas 

We derived the oxygen abundance using two calibrations present 
in the literature: the Kobulnicky et al. d 1 9991 ) calibration, valid 
for HII regions, and the Storchi-Bergmann et al. (119981 ) calibra- 
tion, derived in the assumption of photoionization by a typical 
AGN continuum (Mathews & Ferland |I987l ). The first approach 
was also adopted by Athey & Bregman (2009) to derive nebu- 
lar metallicities in 7 early-type galaxies, and is probably reason- 
able in the case of photoionization by PAGB stars. The second 
calibration is valid in the case of AGN excitation, and may be 
reasonable if LINERs are low-accretion rate AGNs. Our result 
is the following: irrespective of the adopted calibration, the neb- 
ular oxygen metallicities are significantly lower than the stellar 
ones. The main caveat is that power-law ionizing spectra tend 
to provide abundances up to ~0.5 dex larger compared to a typ- 
ical AGN continuum. A calibration in the case of shock heat- 
ing is not available in the literature. However, we can see that 
shock models indicate nebular metallicities not higher than so- 
lar in our sample, in agreement with the results obtained with 
the Kobulnicky et al. (1 1 9991 ) and Storchi-Bergmann et al. d 1 9981) 
calibrations. 

Our result has an important implication for the origin of the 
gas, suggesting that it may have at least in part an external ori- 
gin (e.g., from a cooling flow or from accretion from a com- 
panion galaxy). This is in agreement with photometric and kine- 
matical studies in ETGs showing in many cases gas/star mis- 
alignment and gas/star angular momentum decoupling (Bertola 
et al. [19921 van Dokkum & Franx ,1995. Caon et al. 2000; Sarzi 
et al. 120061 ). Indeed, kinematical and morphological peculiari- 
ties are present in 50% of our sample, suggesting that accretion 
is frequent. Accretion of fresh gas can feed the central supermas- 
sive BH, and activate the AGN emission. This is consistent with 
the fact that a major fraction of our galaxies studied in the X-ray 
and/or radio show nuclear compact sources. 

Finally, we can not exclude that the low oxygen abundance 
is unrelated to the gas origin, and reflects the fact that oxygen 
does not vary on lockstep with the other c-elements. In a sur- 
vey of 27 Milky Way bulge giants, Fulbright, McWilliam, & 
Rich (120061 ) found that, although Mg appears to be enhanced 
at all [Fe/H], [O/Fe] declines with increasing [Fe/H] and is solar 
or mildly sub-solar at [Fe/H] > 0. Also, X -ray studies have de- 
rived subsolar [O/Fe] values, but supersolar [Mg/Fe] values, in 
the hot gas of elliptical galaxies (e.g., Humphrey & Buote l2006l 
Ji et al. 2009). This may be the evidence of an overestimate of 
O yield by SNII models, which do not consider significant mass 
loss at the late stage of massive progenitor stars (Ji et al. ,2009) . 
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Table 2. Observed emission line fluxes in units of 10 '^erg s ' cm - arcsec 



Galaxy 


Regi 


E(B- 


V) 


[O 11]^ 


[Ne III] 


H/3 


[O 111] 


[O 111] 


[O 1] 


[Nil] 


Ho- 


[Nil] 


[SII] 


[SII] 












xiin 


^3869 


i4861 


^4959 


i5007 


.16300 


i6548 


i6563 


/i6584 


A61\l 


i6731 


NGC128 


1 


0.26 


± 


0.29 


6.11 ± 0.55 


0.40 ± 0.55 


1.76 ±0.50 


1.35 ± 0.20 


2.71 ± 0.25 


0.99 ± 0.06 


2.96 ±0.16 


7.13 ± 0.34 


8.72 ±0.10 


2.70 ± 0.23 


1.51 ± 0.18 


NGC128 


2 


0.03 


± 


0.28 


5.87 ± 0.49 


0.07 ± 0.44 


1.68 ± 0.45 


0.85 ± 0.15 


1.96 ±0.19 


0.51 ± 0.09 


2.29 ±0.17 


5.26 ± 0.30 


6.02 ± 0.04 


2.16 ±0.19 


1.64 ±0.18 


NGC128 


3 


0.15 


± 


0.41 


2.96 ± 0.34 


-0.06 ±0.31 


0.74 ± 0.30 


0.48 ± 0.08 


1.02 ±0.11 


-0.67 ± 0.06 


1.09 ±0.11 


2.68 ± 0.20 


2.99 ± 0.04 


1.07 ± 0.12 


0.98 ± 0.12 


NGC128 


4 


0.03 


± 


0.46 


0.62 ±0.15 


-0.09 ± 0.15 


0.39 ± 0.17 


0.22 ± 0.05 


0.35 ± 0.06 


0.57 ± 0.03 


0.46 ± 0.05 


0.76 ±0.10 


-0.12 ±0.01 


0.36 ± 0.07 


0.39 ± 0.07 



NGC777 5 
NGC777"^ 

NGC7775 
NGC7775 



1 


2.81 


± 


3.53 


0.89 ±0.16 


0.13 ±0.23 


0.03 ± 0.12 


-0.03 ± 0.00 


1.14 ±0.07 


0.81 ± 0.02 


1.13 ± 0.12 


1.71 ± 0.05 


1.51 ± 0.02 


7.44 ±0.13 




2 


0.52 


± 


0.71 


0.46 ±0.10 


-0.21 ± 0.15 


0.16 ±0.12 


-0.09 ± 0.05 


0.67 ± 0.08 


0.15 ±0.03 


0.65 ± 0.09 


0.85 ± 0.04 


0.58 ± 0.01 


3.82 ± 0.09 




3 


0.05 


± 


0.24 


0.20 ± 0.04 


-0.14 ±0.04 


0.14 ±0.03 


-0.00 ± 0.00 


0.23 ± 0.02 


0.16 ±0.02 


0.09 ± 0.01 


0.35 ± 0.01 


0.17 ±0.00 


1.74 ±0.01 




4 


0.05 


± 


0.94 


0.13 ± 0.04 


-0.01 ± 0.05 


-0.04 ± 0.04 


0.02 ± 0.01 


0.06 ± 0.02 


0.13 ± 0.01 


0.04 ± 0.01 


0.07 ± 0.02 


0.00 ± 0.00 


0.43 ± 0.01 





NGC1052 


1 


0.22 


± 


0.04 


NGC1052 


2 


0.38 


± 


0.09 


NGC1052 


3 


0.24 


± 


0.20 


NGC1052 


4 


0.03 


± 


0.45 



NGC1209 


1 


0.04 ± 


0.23 


NGC1209 


2 


0.04 ± 


0.24 


NGC1209 


3 


0.04 ± 


0.25 


NGC1209 


4 


0.04 ± 


0.52 



88.41 ± 0.87 
30.98 ± 0.76 
6.71 ±0.32 
0.60 ±0.14 

12.74 ± 0.96 
7.80 ±0.71 
2.83 ± 0.47 
0.14 ±0.18 



10.79 ± 1.11 
2.71 ± 0.54 
0.53 ± 0.28 
0.04 ± 0.19 

-0.39 ± 0.96 
-0.57 ± 0.84 
-0.63 ± 0.63 
-0.13 ± 0.26 



29.21 : 
7.68 d 



: 1.18 
0.65 



15.98 ±0.55 
8.54 ± 0.56 



73.78 ± 0.64 
19.70 ± 0.39 



53.41 ± 0.45 
13.29 ±0.20 



49.57 : 
15.98 : 



0.94 
0.42 



112.80: 

34.57 ± 



:0.94 

0.53 



1.37 


± 


0.27 


0.74 


± 


0.09 


2.52 


± 


0.13 


1.59 


± 


0.07 


1.98 


± 


0.13 


5.38 


± 


0.18 


0.29 


± 


0.13 


0.07 


± 


0.06 


0.31 


± 


0.08 


0.00 


± 


0.01 


0.34 


± 


0.04 


0.71 


± 


0.08 


3.74 


± 


0.85 


2.04 


± 


0.30 


4.82 


± 


0.37 


1.10 


± 


0.17 


5.46 


± 


0.28 


9.39 


± 


0.46 


3.10 


± 


0.71 


1.72 


± 


0.22 


3.52 


± 


0.29 


0.91 


± 


0.13 


3.43 


± 


0.19 


6.34 


± 


0.36 


1.79 


± 


0.43 


0.80 


± 


0.22 


1.21 


± 


0.28 


0.53 


± 


0.04 


1.08 


± 


0.11 


2.63 


± 


0.21 


0.35 


± 


0.17 


0.31 


± 


0.10 


0.31 


± 


0.18 


0.13 


± 


0.03 


0.26 


± 


0.04 


0.42 


± 


0.09 



NGC1297 


1 


0.52 


± 


0.52 


3.74 ± 


0.43 


0.16 ± 0.51 


0.67 


± 


0.34 


0.64 


± 


0.19 


1.35 


± 


0.24 


0.64 


± 


0.06 


1.38 


± 


0.10 


3.49 


± 


0.23 


NGC1297 


2 


0.06 


± 


0.44 


1.30 ± 


0.17 


-0.01 ± 0.25 


0.40 


± 


0.17 


0.34 


± 


0.10 


0.74 


± 


0.15 


0.41 


± 


0.04 


0.59 


± 


0.04 


1.33 


± 


0.10 


NGC1297 


3 


0.03 


± 


0.50 


0.42 ± 


0.08 


0.07 ±0.10 


0.18 


± 


0.09 


0.19 


± 


0.05 


0.19 


± 


0.05 


0.04 


± 


0.02 


0.19 


± 


0.02 


0.36 


± 


0.05 


NGC1297 


4 


1.00 


± 


1.83 


0.08 ± 


0.02 


0.03 ± 0.04 


0.03 


± 


0.05 


0.03 


± 


0.02 


0.07 


± 


0.02 


0.09 


± 


0.01 


0.14 


± 


0.01 


0.22 


± 


0.02 



109.76 ± 0.55 
32.99 ±0.18 
4.67 ± 0.04 
0.35 ± 0.01 

14.21 ± 0.09 
9.02 ±0.11 

3.07 ± 0.05 
0.43 ± 0.02 

3.04 ± 0.04 

1.08 ±0.02 
0.24 ± 0.01 
-0.03 ± 0.00 



67.68 ± 1.16 
23.96 ± 0.44 
3.62 ±0.15 
0.25 ± 0.05 

4.69 ± 0.36 
2.41 ± 0.25 
0.73 ±0.14 
-0.04 ± 0.05 

2.07 ± 0.15 
0.77 ± 0.06 
0.14 ±0.03 
0.18 ± 0.01 



60.98 ± 
17.44 ± 
2.62 ± 
0.21 ± 

3.94 
2.36 
0.75 ± 
-0.00 ± 

1.91 
0.64 
0.11 
0.04 



NGC1380 


1 


0.25 


± 


0.35 


3.20 


± 


0.84 


-0.58 ± 


0.86 


2.31 


± 


0.80 


1.23 


± 0.25 


3.41 


± 


0.39 


1.85 ± 0.16 


5.32 


± 


0.31 


9.21 


± 


0.47 


12.32 


±0.11 


3.49 


± 


0.37 


3.55 


± 


0.40 


NGC1380 


2 


0.51 


± 


0.47 


2.86 


± 


0.70 


0.00 ± 


0.73 


1.45 


± 


0.68 


0.76 


± 0.20 


2.13 


± 


0.29 


0.94 ±0.15 


3.09 


± 


0.22 


7.44 


± 


0.37 


8.29 


± 0.10 


2.73 


± 


0.29 


1.98 


± 


0.25 


NGC1380 


3 


0.44 


± 


0.49 


1.52 


± 


0.45 


-0.10 ± 


0.50 


0.91 


± 


0.44 


-0.02 


± 0.06 


1.02 


± 


0.22 


0.59 ±0.11 


1.69 


± 


0.14 


4.38 


± 


0.26 


4.66 


± 0.05 


1.80 


± 


0.18 


1.16 


± 


0.16 


NGC1380 


4 


0.02 


± 


0.44 


0.47 


± 


0.23 


0.01 ± 


0.32 


0.61 


± 


0.26 


0.03 


± 0.10 


0.62 


± 


0.15 


0.01 ± 0.02 


0.93 


± 


0.07 


1.64 


± 


0.13 


1.63 


± 0.03 


0.55 


± 


0.09 


0.26 


± 


0.07 


NGC1453 


1 


0.11 


± 


0.19 


10.06 


1 ± 


: 0.51 


0.08 ± 


0.42 


2.72 


± 


0.52 


1.72 


± 0.24 


3.93 


± 


0.27 


1.42 ± 0.12 


4.37 


± 


0.23 


9.21 


± 


0.31 


12.56 


± 0.05 


6.28 


± 


0.29 


4.89 


± 


0.26 


NGC1453 


2 


0.11 


± 


0.21 


5.37 


± 


0.39 


-0.14 ± 


0.38 


1.88 


± 


0.38 


1.35 


± 0.13 


2.59 


± 


0.17 


0.59 ± 0.07 


2.45 


± 


0.14 


4.96 


± 


0.23 


6.10 


± 0.04 


2.65 


± 


0.18 


2.40 


± 


0.20 


NGC1453 


3 


0.11 


± 


0.23 


1.39 


± 


0.16 


-0.17 ± 


0.21 


0.75 


± 


0.17 


0.54 


± 0.06 


0.98 


± 


0.09 


0.09 ± 0.02 


0.58 


± 


0.06 


1.34 


± 


0.10 


1.43 


±0.02 


0.40 


± 


0.06 


0.59 


± 


0.07 


NGC1453 


4 


0.11 


± 


0.29 


0.43 


± 


0.06 


-0.03 ± 


0.09 


0.21 


± 


0.06 


0.09 


± 0.03 


0.36 


± 


0.07 


-0.00 ± 0.02 


0.15 


± 


0.02 


0.44 


± 


0.03 


0.40 


± 0.01 


0.14 


± 


0.02 


0.08 


± 


0.02 


NGC15215 


1 


1.22 


± 


1.51 


0.80 


± 


0.40 


0.04 ± 


0.63 


0.24 


± 


0.36 


0.51 


± 0.12 


0.49 


± 


0.23 


-0.00 ± 0.03 


1.16 


± 


0.14 


2.54 


± 


0.25 


3.27 


± 0.06 


0.56 


± 


0.22 


0.46 


± 


0.19 


NGC15215 


2 


0.30 


± 


1.20 


0.54 


± 


0.37 


0.05 ± 


0.59 


0.30 


± 


0.36 


0.07 


± 0.18 


0.17 


± 


0.28 


-0.01 ± 0.02 


0.33 


± 


0.09 


1.26 


± 


0.19 


1.40 


± 0.05 


0.48 


± 


0.09 


0.28 


± 


0.07 


NGC15215 


3 


0.04 


± 


8.66 


0.32 


± 


0.17 


-0.29 ± 


0.27 


-0.02 


± 


0.18 


-0.01 


±0.10 


0.14 


± 


0.15 


0.00 ± 0.01 


0.42 


± 


0.04 


0.29 


± 


0.09 


0.52 


± 0.03 


0.10 


± 


0.04 


0.09 


± 


0.04 


NGC15215 


4 


0.04 


± 


0.77 


0.04 


± 


0.05 


-0.06 ± 


0.08 


0.08 


± 


0.06 


0.08 


± 0.04 


0.22 


± 


0.06 


0.03 ± 0.01 


-0.01 


± 


0.01 


0.14 


± 


0.04 


0.10 


± 0.01 


0.01 


± 


0.01 


0.01 


± 


0.01 



: 0.99 
: 0.32 

0.12 > 

0.03 2. 
a- 
B. 

0.31 

0.22 ^ 

0.11 r 

±0.04 2; 

tO.ll a- 
t 0.04 CD 
1 0.01 B. 
tO.Ol t 
Td 

CO 



3- 
o' 

n' 
ct 
o. 



NGC1533 


1 


0.01 


± 


0.24 


9.45 


± 


1.24 


0.43 ± 0.76 


4.77 ± 1.13 


2.71 


±0.40 


7.99 


± 


0.47 


3.88 ± 0.19 


6.02 


± 0.32 


12.56 ± 0.63 


13.90 ± 0.09 


4.66 


± 


0.42 


3.82 


± 


0.38 


NGC1533 


2 


0.01 


± 


0.23 


4.99 


± 


0.70 


-0.42 ± 0.80 


2.94 ± 0.65 


1.31 


±0.23 


3.54 


± 


0.30 


1.23 ± 0.13 


2.06 


± 0.17 


4.71 ±0.36 


4.70 ± 0.08 


1.62 


± 


0.24 


1.34 


± 


0.20 


NGC1533 


3 


0.01 


± 


0.30 


2.54 


± 


0.32 


0.01 ± 0.39 


1.01 ±0.28 


0.83 


± 0.12 


1.90 


± 


0.18 


0.29 ± 0.06 


0.56 


± 0.07 


1.72 ±0.15 


1.37 ±0.03 


0.88 


± 


0.10 


0.65 


± 


0.08 


NGC1533 


4 


0.01 


± 


0.45 


0.62 


± 


0.11 


-0.01 ± 0.14 


0.18 ± 0.08 


0.17 


± 0.05 


0.30 


± 


0.04 


0.01 ± 0.01 


0.16 


± 0.02 


0.50 ± 0.05 


0.39 ± 0.01 


0.33 


± 


0.03 


0.17 


± 


0.02 



3 Reg 1 = r <r,/16; Reg 2 = r,/16 < r <r,/8; Reg 3 = r,/8< r <r,/4; Reg 4 

Negative emission values are consistent with no emission. 
^ The continuum was fitted with 2 SSPs 



= r^/4 < r <rJ2. 



Table 3. Observed emission line equivalent widths. 



Galaxy 




[oiiii 


[Ne III] 


Ufi 


[O III] 


[O III] 


[O I] 


[Nil] 


Hff 


[Nil] 


[SII] 


[SII] 






i3727 


^3869 


i4861 


i4959 


A5007 


^6300 


^6548 


i6563 


^6584 


i6717 


/i6731 



NGC128 


I 


7.25 ± 0.59 


0.59 


± 


0.81 


0.57 


± 0.16 


0.43 


± 0.06 


0.87 


± 


0.08 


0.27 


± 0.02 


0.74 


± 


0.04 


1.78 


± 


0.09 


2.18 ± 0.03 


0.73 


± 


0.06 


0.41 


± 


0.05 


NGC128 


2 


7.13 ±0.59 


0.13 


± 


0.75 


0.64 


± 0.17 


0.32 


± 0.05 


0.73 


± 


0.07 


0.17 


± 0.03 


0.69 


± 


0.05 


1.58 


± 


0.09 


1.81 ± 0.01 


0.71 


± 


0.06 


0.54 


± 


0.06 


NGC128 


3 


5.00 ± 0.58 


-0.15 


± 


0.74 


0.40 


± 0.16 


0.25 


± 0.04 


0.54 


± 


0.06 


-0.31 


±0.03 


0.48 


± 


0.05 


1.18 


± 


0.09 


1.32 ±0.02 


0.51 


± 


0.06 


0.47 


± 


0.06 


NGC128 


4 


1.92 ± 0.48 


-0.41 


± 


0.68 


0.39 


± 0.17 


0.22 


± 0.05 


0.34 


± 


0.06 


0.51 


±0.03 


0.38 


± 


0.04 


0.62 


± 


0.08 


-0.10 ±0.01 


0.32 


± 


0.06 


0.34 


± 


0.06 


NGC777'* 


1 


1.36 ± 0.24 


0.27 


± 


0.50 


0.01 


± 0.04 


-0.01 


± 0.00 


0.42 


± 


0.03 


0.27 


± 0.01 


0.36 


± 


0.04 


0.54 


± 


0.02 


0.47 ± 0.01 


2.53 


± 


0.05 










2 


1.16 ± 0.25 


-0.71 


+ 


0.50 


0.11 


± 0.08 


-0.06 


± 0.03 


0.47 




0.06 


0.09 


± 0.02 


0.37 




0.05 


0.49 




0.02 


0.33 ± 0.01 


2.38 




0.06 








NGC777^ 


3 


1 10 + 90 


-1 04 




0.28 


23 


+ 05 


-0 00 


+ 00 


0.39 




0.04 


24 


+ 03 


12 




02 


0.49 




02 


94 + 01 


2.64 




02 








NGC777'* 


4 


2.34 ± 0.74 


-0.31 


± 


1.10 


-0.22 


± 0.20 


0.09 


±0.05 


0.33 


± 


0.12 


0.66 


±0.03 


0.17 


± 


0.04 


0.33 


± 


0.08 


0.02 ± 0.01 
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^ Reg 1 = r <r,/16; Reg 2 = r,/16 < r <r,/8; Reg 3 = r,/8< r <r,/4; Reg 4 
^ Negative emission values are consistent with no emission. 
^ The continuum was fitted with 2 SSPs 



= r^/4 < r <rJ2. 
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Appendix A: Short description of morphological 
and kinematic characteristics and peculiarities 
of the galaxy sample (Paper l+Paper II) 

We collect in Table Al the morphological and kinematic charac- 
teristics, as well as the peculiarities, of the complete R05-(-A06 
sample. The morphological and kinematic properties in Table Al 
are derived from the literature, and were in part already de- 
scribed in the Appendix (provided as on-line material) of Paper I 
and Paper II. 

We list in Cols. 3 and 4 the galaxy ellipticity, e, and 
the {V I cr) scaled- The latter is computed from slit spectroscopy 
measurements through the relation provided by Cappellari et 
al. (120071) . (y/o-)scaied =0.57x(V/cr)^,„. It should be a reason- 
able approximation of the (V/cr)^, computed on a luminosity 
weighted spectrum within r^. For the eight galaxies in common 
with the SAURON sample, i.e. NGC 2974, NGC 3489, NGC 4374, 
NGC 4552, NGC 5813, NGC 5831, NGC 5846 and NGC 7332, 
the e and V/cr values are those derived through integral field 
spectroscopy. Cappellari et al. ( I2007I I noticed that there is a good 
agreement between e values derived in the literature with differ- 
ent methods. 

Emsellem et al. (120071) defined a new parameter, 
Ar={R\V\}/{R VV^ + 0-2), where R is the galacto-centric 
distance, and V and cr are luminosity weighted averages over 
the two-dimensional kinematic field. This parameter is a proxy 
to quantify the observed projected stellar angular momentum 
per unit mass, and can be used to classify ETGs into slow or 
fast rotators. In the e vs. (yicr)e plane, Cappellari et al. ( I2007I I 
identified a region (0.0< e <0.3, 0.0< (V/cr), <0.13) where 
slow- rotators are generally found. Lacking the 2D V and cr 
information for most of our galaxies, we used {'V/(T)scaied to 
identify in the e vs. {'V/cr)g plane the fast and slow rotators. The 
derived classification is provided in Col. 5 of Table Al (F=fast 
rotator; S=slow rotator). We obtain that x6E% (36/53) of our 
sample is fast rotator (to be compared with 75% fast rotators 
in the SAURON sample). Less than 1/3 of our galaxies are slow 
rotators (vs 25% in the SAURON sample). 

In Col. 6 and 7 of Table Al we provide the kinematic and 
morphological peculiarities, respectively. A large fraction, ap- 
proximately 50% of our ETGs, displays kinematic peculiarities, 
like counter-rotation, or peculiar faint features, like shells, tails, 
or decoupled gas distributions with respect to the stars. 
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Table Al. kinematic and morphological overview 



ident 


RSA 


e 


{WI(T)scaled 


Rot. 


Gas vs. stars kinematic 


Gas vs. stars morphological 










class 


peculiarities 


peculiarities 


NGC 128 


S02(8) pec 


0.67 


0.44 


F 


CR g-s(l) 


g- d and g-maj tai23° (5) 


NGC 777 


El 


0.21 


0.06 


S 






NGC 1052 


E3/S0 


0.28 


0.27 


F 


CR g-g (1,5) 




NGC 1209 


E6 


0.52 


0.50 


F 




X- like struct. ; NW linear feature (2) 


NGC 1297 


S02/3(0) 


0.13 










NGC 1366 


E7/S01(7) 


0.56 










NGC 1380 


S03(7)/Sa 


0.41 


0.59 


F 






NGC 1389 


S01(5)/SB01 


0.37 










NGC 1407 


E0/S01(0) 


0.07 


0.11 


S 


rotat. min. axis (5) 




NGC 1426 


E4 


0.34 


0.40 


F 






NGC 1453 


EO 


0.17 








g-d and g-maj t (5) 


NGC 1521 


E3 


0.35 


0.54 


F 






NGC 1533 


SB02(2)/SBa 


0.19 










NGC 1553 


S01/2(5)pec 


0.38 


0.68 


F 




shells (3) 


NGC 1947 


S03(0) pec 


0.11 


0.23 


F 


rotat. min. axis (5) 


dust-lane min. axis (5) 


NGC 2749 


E3 


0.07 


0.21 


F 


rotat. min. axis (5) 




NGC 2911 


SOp or S03(2) 


0.32 










NGC 2962 


RSB02/Sa 


0.37 


0.51 


F 






NGC 2974 


E4 


0.37 


0.70 


F 




shells (2); g-d and g-maj tsj20° (5) 


NGC 3136 


E4 


0.24 


0.67 


F 


CR s-s (5) 




NGC 3258 


El 


0.13 


0.08 


S 


CR g-s (1) 




NGC 3268 


E2 


0.24 


0.53 


F 






NGC 3489 


S03/Sa 


0.37 


0.29 


F 


gas rot. min. axis (5) 




NGC 3557 


E3 


0.21 


0.30 


F 




SW fan; asym. outer isophotes (2) 


NGC 3607 


S03(3) 


0.11 


0.28 


F 






NGC 3818 


E5 


0.36 


0.27 


F 






NGC 3962 


El 


0.22 


0.23 


F 




gas disk+outer arc-like struct. (5) 


NGC 4374 


El 


0.13 


0.03 


S 






NGC 4552 


SOl(O) 


0.06 


0.05 


S 




shells (3) 


NGC 4636 


E0/S01(6) 


0.24 


0.07 


S 


gas irr. motion (5) 




NGC 4696 


(E3) 


0.34 


0.18 


F 




Faint outer shells (2) 


NGC 4697 


E6 


0.32 


0.39 


F 




Non spherical inner isophotes (2) 


NGC 5011 


E2 


0.15 


0.11 


S 






NGC 5044 


EO 


0.11 


0.13 


F 


CR s-s; gas irr.motion (5) 


gas lil. shape (5) 


NGC 5077 


SO 1/2(4) 


0.15 


0.06 


S 


CR s-s (5); g-d and g-maj tai90° (5) 




NGC 5090 


E2 


0.15 


0.17 


F 






NGC 5193 


SO 1(0) 


0.07 










NGC 5266 


S03(5) pec 


0.31 


0.12 


S 


stellar rot. min. axis (5) 


dust-lane and gas along min. axis (5) 


NGC 5328 


E4 


0.31 










NGC 5363 


[S03(5)] 


0.34 


0.39 


F 




dust -lane along min. axis (5) 


NGC 5638 


El 


0.11 


0.21 


F 






NGC 5812 


EO 


0.13 


0.11 


S 




Tidal tail (2) 


NGC 5813 


El 


0.15 


0.14 


S 


CR s-s; gas irr. motion (5) 


gas til. shape (5) 


NGC 5831 


E4 


0.15 


0.08 


s 


CR s-s (6,7) 




NGC 5846 


SO 1(0) 


0.07 


0.03 


s 


gas iiT. motion (5) 


Faint outer shells (2) 


NGC 5898 


S02/3(0) 


0.07 


0.21 


F 


CR g-s(l) 


Three spiral arm-like tidal tails (2) 


NGC 6721 


El 


0.15 


0.26 


F 






NGC 6758 


E2 (merger) 


0.22 


0.12 


S 


CR g-s (5) 




NGC 6776 


El pec 


0.17 


0.24 


F 




shell or loops NE (3,4); gas til. shape (5) 


NGC 6868 


E3/S02/3(3) 


0.19 


0.14 


F 


CR g-s (5) 




NGC 6875 


SO/a(merger) 


0.41 










NGC 6876 


E3 


0.13 










NGC 6958 


R?S01(3) 


0.15 


0.20 


F 




Shells (2,3) 


NGC 7007 


S02/3/a 


0.42 






CR g-s (1) 


g-d and g- maj t»30° (5) 


NGC 7079 


SBa 


0.32 


0.96 


F 


CR g-s + s-s (1) 




NGC 7097 


E4 


0.29 


0.10 


S 


CR s-s (1) 




NGC 7135 


SOI pec 


0.31 


0.20 


F 




jet-like debris and shells (3) 


NGC 7192 


S02(0) 


0.15 


0.04 


S 


CR s-s (5) 


Shell (2) 


NGC 7332 


S02/3(8) 


0.42 


0.32 


F 


CR g-g (1) 




NGC 7377 


S02/3/Sa pec 


0.19 










IC 1459 


E4 


0.28 


0.10 


S 


CR g-g (1) 


shells (2) 


IC 2006 


El 


0.15 


0.16 


F 


CR g-s(l) 




IC 3370 


E2 pec 


0.21 


0.28 


F 




X-like struct.; broad N fan (2); polar ring ? (5) 


IC 4296 


EO 


0.17 


0.13 


F 


CR s-s (5) 




IC 5063 


S03(3)pec/Sa 


0.28 


0.26 


F 







Cols. 4 and 5 list respectively (V/(T)jra/(,rf=0.57(V/(T)j/,f, and the rotation class (F=fast rotator; S=slow rotator). Legenda: CR g-s: counter rotation gas vs. stars; 
CR s-s: counter rotation stars vs. stars; CR g-g: counter rotation gas vs. gas; stars rotat. min. axis: stars rotate along the galaxy minor axis; g-d and g-maj t: gas 
disk and galaxy major axis ai'e tilted by the reported angle, if provided in the literature. References: (1) Corsini & Bertola 1998, (2) Tal et al. 2009, (3) Malin & 
Carter ! 1 983 i (4) Pierfederici & Rampazzo i2004i (5) the description of the kinematic and morphological peculiarities of the galaxies and full references are reported 
in the on-line notes of Paper I and Paper II; (6) Davies et al. ll983l (7) Emsellem et al. l2004l 
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Appendix B: Emission line ratios for the total 
sample in annuli of increasing galacto-centric 
distance. 
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Fig. B.l. log[Nn]/Ha ratio for the sample in the annuh r <re/16, r^/ie < r <re/8, re/8< r <re/4, and re/4 < r <re/2. 
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Fig. B.2. log[Nn]/Ha ratio for the sample in the annuh r <re/16, r^/ie < r <re/8, re/8< r <re/4, and re/4 < r <re/2. 
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Fig. B.3. log[OIII]/Hy8 ratio for the sample in the atmuli r <re/16, r^/ie < r <re/8, re/8< r <re/4, and r^M < r <re/2. 
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Fig. B.4. log[OIII]/Hy8 ratio for the sample in the atmuli r <re/16, r^/ie < r <re/8, re/8< r <re/4, and r^M < r <re/2. 



